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^ We present a method for calculating transfer matrices for the g-state Potts model 
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partition functions Z(G,q,v), for arbitrary q and temperature variable v, on strip 
graphs G of the square (sq), triangular (tri), and honeycomb (he) lattices of width 
L y vertices and of arbitrarily great length L x vertices, subject to toroidal and Klein 
bottle boundary conditions. For the toroidal case we express the partition function 
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as Z(A,L y x L x ,q,v) = J2d^o^j b j i^Z,A,L v ,d,j) m , where A denotes lattice type, 



bj are specified polynomials of degree d in q, Xz,A,L y ,d,j are eigenvalues of the 

in ; 

transfer matrix Tz : \,L y ,d in the degree-d subspace, and m = L x (L x /2) for A = 
sq, tri (he), respectively. An analogous formula is given for Klein bottle strips. We 
exhibit a method for calculating Tz a,l d f° r arbitrary L y . In particular, we find 
some very simple formulas for the determinant det(Tz t A,L y ,d), and trace Tr(Tz,A,L y )- 



Corresponding results are given for the equivalent Tutte polynomials for these lattice 
strips and illustrative examples are included. 
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I. INTRODUCTION 



The g-state Potts model has served as a valuable model for the study of phase transitions 
and critical phenomena On a lattice, or, more generally, on a (connected) graph G, 

at temperature T, this model is defined by the partition function 



with the (zero-field) Hamiltonian 



H=-jJ2Sa i a j (1-2) 



where crj = 1, ...,q are the spin variables on each vertex (site) i 6 G; (3 — (/c^T) -1 ; and (ij) 
denotes pairs of adjacent vertices. The graph G = G(V, E) is defined by its vertex set V and 
its edge set E\ we denote the number of vertices of G as n = n(G) = \V\ and the number of 
edges of G as e(G) = \E\. We use the notation 



K = pJ , v = e K - 1 (1.3) 

so that the physical ranges are v > for the Potts ferromagnet, and — 1 < v < for Potts 
antiferromagnet, corresponding to < T < oo. One defines the (reduced) free energy per 
site / = —f3F, where F is the actual free energy, via f({G}, q, v) = lim^oo \n[Z(G, q, v ) 1//n ], 
where we use the symbol {G} to denote the formal limit lim^oo G for a given family of 
graphs. In the present context, this n — > oo limit corresponds to the limit of infinite length 
for a strip graph of fixed width and some prescribed boundary conditions. 

In this paper we shall present transfer matrices for the g-state Potts model partition 
functions Z(G, q, v), for arbitrary q and temperature variable v, on toroidal and Klein bottle 
strip graphs G of the square, triangular, and honeycomb lattices of width L y vertices and of 
arbitrarily great length L x vertices. We label the lattice type as A and abbreviate the three 
respective types as sq, tri, and he. Each strip involves a longitudinal repetition of m copies 
of a particular subgraph. For the square-lattice strips, this is a column of squares. It is 
convenient to represent the strip of the triangular lattice as obtained from the corresponding 
strip of the square lattice with additional diagonal edges connecting, say, the upper-left to 
lower-right vertices in each square. In both these cases, the length is L x = m vertices. We 
represent the strip of the honeycomb lattice in the form of bricks oriented horizontally. In 
this case, since there are two vertices in 1-1 correspondence with each horizontal side of a 
brick, L x = 2m vertices. Summarizing for all of three lattices, the relation between the 
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number of vertices and the number of repeated copies is 



, m if A = sq or tri 
L x = { . (1.4) 

2m if A = he 



For the toroidal case the partition function has the general form 

Ly 

Z(A,L y x L x ,tor.,q,v) =£^if(A W3 ) m (1.5) 

d=0 j 

where Xz,A,L y ,d,j, are eigenvalues of the transfer matrices in the degree-d subspace T Zt \,L y ,d 
and they are independent of the length of the strip length L x . The size of the transfer 
matrices T Z) A,L y ,d is equal to d\ n ztor (L yi d) that will be discussed in the following section. 
Here because the dimensions of Tz,h,L y ,d are the same for all three lattices A = sq, tri, he, we 
omit A in the notation. In contrast to cyclic and Mobius strips, there can be more than one 
coefficient, denoted as bj, for each degree d. The fe^'s are polynomials of degree d in q and 
play a role analogous to multiplicities of eigenvalues Xz,A,L y ,d,j, although this identification is 
formal, since bj may be zero or negative for the small physical values of q. We exhibit our 
method for calculating Tz t A,L y ,d f° r arbitrary L y , and we shall construct an analogous formula 
for Klein bottle strips. Explicit results for arbitrary L y are given for (i) Tz t A,L y ,d with d = L y , 
(ii) the determinant det(Tz t A,L y ,d), and (iii) the trace Tr(Tz,A,L v )- The results in (i) and (iii) 
apply for A = sq, tri, he, while the results in (ii) apply for A = sq, he. Corresponding results 
are given for the equivalent Tutte polynomials for these lattice strips and illustrative examples 
are included. We have calculated the transfer matrices up to widths L y = 4 for the square 
and honeycomb lattices and up to L y = 3 for the triangular lattice. Since the dimensions 
of these matrices increase rapidly with strip width (e.g., dim{T z ,A,L y ,d) = 14,35,56,48 for 
L y = 4 and < d < 3), it is not feasible to present many of the explicit results here; instead, 
we concentrate on general methods and results that hold for arbitrary L y . 

Various special cases of the Potts model partition function Z{G, q, v) are of interest. For 
example, if one considers the case of antiferromagnetic spin-spin coupling, J < and takes 
the temperature to zero, so that K = — oo and v = —1, then 

Z(G,q,-l)=P(G,q) (1.6) 

where P(G, q) is the chromatic polynomial (in q) expressing the number of ways of coloring 
the vertices of the graph G with q colors such that no two adjacent vertices have the same 
color (5, 0|. The minimum number of colors necessary for such a coloring of G is called the 
chromatic number, x(G)- The (7- state Potts antiferromagnet (AF) 0, 0] exhibits nonzero 
ground state entropy, Sq > (without frustration) for sufficiently large q on a given lattice A 
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or, more generally, on a graph G = (V, E). This is equivalent to a ground state degeneracy 
per site W > 1, since S = k B \nW. Thus 

W({G},q)= hmP(G,g) 1 /". (1.7) 

n —>oo 

For a strip graph G s of a lattice A with given boundary conditions, following our earlier 
notation |3j we denote the sum of coefficients (generalized multiplicities) cz,gj of the Xz,G,j 
as 

C z ,g = Yl Cz &o (!-8) 
i=i 

while for the chromatic polynomial as 

C P ,G = Yl C ^,j ■ (1-9) 
J'=l 

These sums are independent of the length m of the strip. General results for these sums will 
be given below for the strips of interest. 

We recall some previous related work. The partition function Z(G,q,v) for the Potts 
model was calculated for arbitrary q and v on strips of the lattice A with toroidal or Klein 
bottle longitudinal boundary conditions was calculated for (i) A = sq, L y — 2, 3 in 0], 
(ii) for L y = 3 on the square-lattice with next-nearest-neighbor spin-spin couplings 
Matrix methods for calculating chromatic polynomials were developed and used in 
12, is| and more recently in [lj, EE [3]. Ref. 17, 18, 3| developed transfer matrix 
methods for both Z(G, q, v) and the special case v = —1 of chromatic polynomials on strips of 
the square and triangular lattices with free longitudinal boundary conditions and used them 
to calculate the latter polynomials for a large variety of widths. These have been termed 
transfer matrices in the Fortuin-Kasteleyn representation (see eq. (11.10)1 below). These 
methods were applied to calculate the full Potts model partition function for strips of the 
square and triangular lattices with free boundary conditions and a number of widths in Refs. 
|20Ll21|. and then were generalized for strips of the square, triangular and honeycomb lattices 
with cyclic and Mobious strips [i^, 23]. Here we extend these transfer matrix methods for 
calculating Potts model partition functions on toroidal and Klein bottle lattice strips and 
present general results for these lattice strips of arbitrary width. Clearly, new exact results 
on the Potts model are of value in their own right. 

Let G' = (V, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex 
set V and an edge set E' C E. Z(G, q, v) can be written as the sum 0, 

Z(G,q,v)= J2^ G '^ E ' 1 (1-10) 

G'CG 
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where k(G') denotes the number of connected components of G' . Since we only consider 
connected graphs G, we have k(G) = 1. The formula f|l . 1(J|) enables one to generalize q from 
Z + to 5Lf for physical ferromagnetic v. More generally, eq. (jl.lOj) allows one to generalize 
both q and v to complex values, as is necessary when studying zeros of the partition function 
in the complex q and v planes. 

The Potts model partition function Z(G,q,v) is equivalent to an object of considerable 
current interest in mathematical graph theory, the Tutte polynomial, T(G,x,y), given by 



31. 32] 



T(G, x,y)= 2 (x - i) k (G>)-HG)( y _ 1)C (G') (LU) 

G'CG 



where c(G') = \E'\ + k(G') — |V| is the number of independent circuits in G 1 . Now let 

x = l + -, y = v + l (1.12) 

v 



so that 



q=(x-l){y-l) . (1.13) 



Then the equivalence between the Potts model partition function and the Tutte polynomial 
for a graph G is 

Z(G,q,v) = (x-l) k ^(y-ir^T(G,x,y) . (1.14) 

Given this equivalence, we can express results either in Potts or Tutte form. We will use 
both, since each has its own particular advantages. The Potts model form, involving the 
variables q and v is convenient for physical applications, since q specifies the number of states 
and determines the universality class of the transition, and v is the temperature variable. 
The Tutte form has the advantage that many expressions are simpler when written in terms 
of the Tutte variables x and y. 

From (|1.14j) . one can write the Tutte polynomial as 

T{A,L y x L x ,tar.,x,y) = ^SE 6 f^Aj') m (1.15) 

d=0 j 

where m is given in terms of L x by eq. (jl.4j) and it is convenient to factor out a factor of 
l/(x — 1). (This factor is always cancelled, since the Tutte polynomial is a polynomial in x 
as well as y.) 

From eq. (Oil it follows that H Q 



\Z,A,Ly,d,j 



•"'••A/.A./.,./, (1-16) 
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and 

T z ,A, Ly ,d = v pLv T T;A>Ly>d (1.17) 

where 

{1 if A = so or tri 
(1.18) 
2 if A = he 

Note that the factor of {x — 1) in eq. (|1.14j) cancels the factor — 1) in eq. ([1.150 . 



II. TRANSFER MATRIX METHOD 

The chromatic polynomials for strips with periodic boundary condition in the longitudinal 
direction were discussed in Refs. |l4l llfiL llfil ] in terms of a compatibility matrix, and the bases 
of the matrix are given by [AT|£], where X is a subset of L y vertices and £ is an injection from 
X to {1, 2, q}. We now use the same bases for the transfer matrix of the full Potts model 
partition function. The degree d = subspace of the transfer matrix (equivalently called the 
"level 0" subspace in jig ) corresponds to the empty set X = 0, and the bases of the transfer 
matrix are all of the possible non-crossing partitions of L y vertices. (For the zero-temperature 
antiferromagnetic Potts model, a non-nearest-neighbor requirement is imposed; this will be 
discussed further in the next section.) The dimension of this matrix is n^ jior (L y , 0) = Cl v , 
the Catalan number ^?|, as in eq. (j2.21[) below. The eigenvalues of the transfer matrix 
for a cylindrical strip are a subset of the eigenvalues of the transfer matrix Tz,k,L y ,d=o in 
this degree d — subspace for the corresponding toroidal strip. The degree d = 1 subspace 
of the transfer matrix is given by all of the possible non-crossing partitions with a color 
assignment to one vertex (with possible connections with other vertices); i.e., X contains 
one vertex, and the multiplicity is q — 1 = . This follows because there are q possible ways 
of making this color assignment, but one of these has to be subtracted, since the effect of 
all the possible color assignments is equivalent to the choice of no specific color assignment, 
which has been taken into account in the level subspace. In this derivation and subsequent 
ones we assume that q is a sufficiently large integer to begin with, so that the multiplicities 
are positive-definite; we then analytically continue them downward to apply in the region of 
small q where the coefficients can be zero or negative. For the next subspace we consider 
all of the non-crossing partitions with two-color assignments to two separated vertices (with 
possible connections with other vertices). Now the multiplicity can be understood by the 
sieve formula of 0, [l5, Since the two assigned colors should be different, there are 
q(q — 1) ways of making these assignments. This includes the q possible color assignments 
for each of the two vertices that have been considered in level 1 and hence these must be 
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subtracted. In doing this, the no-color assignment was subtracted twice, and one of these 
has to be added back. Therefore, the multiplicity is 

q(q-l) ~2q + l = q 2 -3q + l = b {2) . (2.1) 

Recall that a slice of cyclic strips is a tree graph, but a slice of toroidal strips is a circuit 
graph. For toroidal strips, these coefficients with < d < 2 are the same as those for 
eyciie strips c« in M . (see also 0) b ecause when the number of vertices is 1 or 2 a 
tree graph is equivalent to a circuit graph modulo multiple edges. For toroidal strips, the 
multiplicity for three-color assignment is 

q(q - l)(q - 2) - 3q(q - 1) + 3q - 1 = q 3 - 6q 2 + 8q - 1 = 6 (3) . (2.2) 

By similar reasoning, the multiplicity for four-color assignment is 

q(q - l)(q 2 - 3q + 3) - Aq(q - I) 2 + Aq(q - 1) + 2q 2 - Aq + 1 

= q 4 - 8q 3 + 20q 2 - 15q + 1 = 6 (4) . (2.3) 
Further, we find 
5(5) = q 5 _ 1Qg 4 + 35q 3 _ 5Qq 2 + 24g - 1 

fe( 6 ) = q 6 - I2q b + 54g 4 - 112g 3 + 105g 2 - 35g + 1 

M 7 ) = q 7 - Uq 6 + 77q 5 - 210g 4 + 294g 3 - 196g 2 + A8q - 1 

feW = ( g 2 - 3g + l)(g 6 - 13g 5 + 64g 4 - 147g 3 + 155g 2 - 60q + 1) 

b (9) = g 9 _ lgq 8 + 135g 7 _ 54gg 6 + 12g7q 5 _ 17g2g 4 + 13g6(? 3 _ 54Qg 2 + gQg _ 1 

6 (io) = g io _ 2Qg 9 + l7Qq 8 _ gQQq 7 + 2275g 6 - 4004g 5 + 4290g 4 - 2640g 3 + 825g 2 - 99g + 1 

(2.4) 

and so forth for higher values of d. We show the calculation for these multiplicities pictorially 
for 2 < d < A in Fig. ^ Let us define 

q = 2 + t + t' 1 = 2 + 2cos£9 = 4cos 2 (£?/2) . (2.5) 
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Then, in terms of these variables, we find that &w can be written as 

6 (0) = 1 (2.6) 
&M =t+l + r 1 = l + 2cos# (2.7) 

and, for d > 2, 

= (-l) d (t + l+t- 1 )+t d + r d = (-l) d (l + 2 cos 9) + 2 cos(d#) (2.8) 
Therefore, we have the recursion relation 
b (d+2) + (_yd+i b (i) = ( f + + (_i)^W) _ + (_i)d-ifoW) f or d > 2 , (2.9) 



that is 

& (d+2) = ( g _ 2)^+i) _ fo(d) + g (_i)^ (1) for d > 2 . 
Now the coefficient 

c (d) in 

|3j| is given by 

d 



i=o 



(2.10) 



(2.11) 



where U n (x) is the Chebyshev polynomial of the second kind. Equivalently, in terms of the 
angle 6 in eq. (J23J) Ql^, 

u) _ sin ((2d +l)f) 



sin (|) 



In general, we find that the relation between £>w and is as follows, 

c {d) = b (d) f or o < d < 2 



(2.12) 



and 



bid ' ] f or odd d 



<2'=1 



c (d) = ^ 6 («0 for 

even d 



d'=2 



6 (<0 = c W_ c (d-D + (_i)d c (D 

= u 2d (y^)-u 2d ^) + (-i) d ( q -i) 



(2.13) 



= / + E(-l) J ^-( 2d . J V j + - 1) forrf>2 (2.14) 

i=i J V J / 

At the value q = 0, few has the property that 

&W( g = 0) = . (2.15) 

We also find the following properties of the few polynomials. For g = 1 and d > 2, few i s 
equal to 2 if d — mod 3 and —1 if d — 1 or d — 2 mod 3. For q = 2 and d > 2, few 
is equal to 3 if d = mod 4 and —1 if d = 1, 2, 3 mod 4. One can derive similar relations 
for other values of g. For q = 3, few takes on the pattern of values —4, 1, —1,4, —1, 1 for 
d = 3,4, 5, 6, 7, 8, and then this pattern repeats for higher integral values of d. For q = 4 and 
d > 2, few = 5 if g? is even and few = — 1 if d is odd. For larger values of q, few does not cycle 
through a fixed pattern of values as d increases, but instead increases. The fixed patterns 
of values of few for integer q from to 4 are related to special sum rules for Potts model 
partition functions on lattice sjrips with toroidal boundary conditions, similar to those that 
we exhibited in Refs. M and 
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O-ffl 

©-© 



? 

C) 

© 



+ 



+ 2 



O 

o 



4 O + 1 



FIG. 1: Multiplicities for 2 < d < 4. 



The partitions for toroidal strips are more involved than those for cyclic strips because 
of the rotational symmetry of the slice of toroidal strips. We list graphically all the possible 
partitions for strips of widths L y = 2 and L y = 3 in Figs. 121 and E3 respectively, where 
white circles are the original L y vertices and each black circle corresponds to a specific 
color assignment. In the following discussion, we will simply use the names white and 
black circles with these meanings understood implicitly. The connections between the black 
circles and white circles should obey the non-crossing restriction. The apparent crossings 
in some partitions in Figs. 121 and H3 do not violate this rule when white circles are rotated 
appropriately, owing to the rotational symmetry for toroidal strips. For example, the crossing 
in the last partition of d = 1 in Fig. El can be removed by moving the bottom white circle 
to the top so that it becomes equivalent to the fourth partition of d — 1, or moving the top 
white circle to the bottom so that it becomes equivalent to the eighth partition of d = 1 in 
the figure. Another example is the crossing in the eleventh partition of d = 2 in Fig. El that 
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d = O 

o 



d= 1 o 



2 



? 
o 




FIG. 3: Partitions for the L y = 3 strip. 



can be removed by the rotational symmetry of the white circles. In addition, it is necessary 
to include all the arrangements of black circles when the transfer matrix is constructed. For 
example, the two partitions of d = 2 in Fig. El correspond to exchanging the black circles, 
i.e. the assignments of the colors. We arrange the partitions of d = 2 in Fig. |3] into pairs 
such that every two partitions are equivalent when the two black circles are exchanged. We 
show the six partitions of d = 3 in Fig. 01 which correspond to all the possible arrangements 
of black circles. 

We denote the partitions VL y ,d for 2 < L y < 4 as follows: 

P 2 ,o = {/;12}, 7> 2) i = {2;!;T2}, V 2 , 2 = {1, 2; i, 2} (2.16) 



P 3i0 = {J; 12; 13; 23; 123} , V 3A = {3; 2; 1; 12, 3; 12; 13; 23; 23, 1; 123; 13, 2} , 
P 3 , 2 = {2,3;2,3;l,3;i ) 3;l,2;i,2;12,3;f2,3;l,2^;i,^ 

P 3 ,3 = {1, 2, 3; i, 2, 3; 1, 2, 3; 1, 2, 3; 1, 2, 3; 1, 2, 3} (2.17) 
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V 4 ,o = {/; 12; 13; 14; 23; 24; 34; 12, 34; 14, 23; 123; 124; 134; 234; 1234} , 

V 4 ,i = {4; 3; 2; 1; 12, 4; 12, 3; T2; 13, 4; 13; Ti; 23, 4; 23; 23, 1; 24; 24, 1; 34; 34, 2; 34, 1; 12, 34; 

34, T2; 23, 14; 123, 4; T23; 124; T34; 234; 234, 1; 1234; 13, 2; 14, 3; 14, 2; 24, 3; 14, 23; 

124, 3; 134, 2} , 

P 4 , 2 = {3, 4; 3, 4; 2, 4; 2, 4; 1,4; i, 4; 2, 3; 2, 3; 1, 3; 1,3; 1, 2; 1, 2; 12, 3, 4; 12, 3, 4; T2, 4; f2, 4; 

12, 3; 12, 3; 13, 4; 13, 4; 13, 2; 13, 2; 14, 3; 14, 3; Ti, 2; 14, 2; 14, 2, 3; 14, 2, 3; 23, 4; 23, 4; 
23, 1, 4; 23, 1, 4; 1, 23; i, 23; 1, 24; i, 24; 24, 3; 24, 3; 2, 34; 2, 34; 1, 34; 1, 34; 34, 1, 2; 
34, i, 2; 12, 34; 12, 34; Ti, 23; 14, 23; 123, 4; 123, 4; T24, 3; 124, 3; 134, 2; 134, 2; 1, 234; 
1,234} , 

P 4 ,3 = {2, 3, 4; 2, 3, 4; 2, 3, 4; 2, 3, 4; 2, 3, 4; 2, 3, 4; 1, 3, 4; 1, 3, 4; 1, 3, 4; i, 3, 4; 1, 3, 4; 1, 3, 4; 

1, 2, 4; 1, 2, 4; 1, 2, 4; 1, 2, 4; 1, 2, 4; 1, 2, 4; 1, 2, 3; i, 2, 3; 1, 2 , 3; 1, 2, 3; 1, 2, 3; 1, 2, 3; 

T2, 3, 4; 12, 3, 4; 12, 3, 4; 12, 3, 4; 12, 3, 4; 1~2, 3, 4; Ti, 2, 3; 14, 2, 3; Ti, 2, 3; 14, 2, 3; 

14, 2, 3; 14, 2, 3; 1, 23, 4; 1, 23, 4; 1, 23, 4; 1, 23, 4; 1, 23, 4; 1, 23, 4; 1, 2, 34; 1, 2, 34; 

1, 2, 34; i, 2, 34; 1,2, 34; 1,2, 34} , 
7\ 4 = {1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 

1, 2, 3, 4; 1,2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 

1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4; 1, 2, 3, 4} (2.18) 

where partitions are separated by a colon, and overline, hat, check and tilde in each partition 
correspond to different color assignments. It follows that the size of the transfer matrix 
Tz,A,L y ,d always has the factor d\. Let us denote the reduced size of the transfer matrix as 
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nz,tor(L y , d), i.e., we neglect the permutation of black circles. As the derivation of does 
not take into account permutations, i.e. a set of d\ eigenvalues should have their coefficients 
summed to be equal to , we find that the sum of all coefficients for a toroidal strip graph, 
which is equal to the dimension of the total transfer matrix, i.e., the number of ways to color 



L y vertices, is 



C z ,A,L y = dim{T ZAiLy ) = n z , tar (L y , d)b {d) = q L y for A = sq, tri, he . (2.19) 



d=0 



By substituting the expression for few in eq. (|2.19|) . we determine the nz,tor(L y , d), as 
follows: 

nz,tor(L y , d) = for d > L y (2.20) 

n z ,tor{L y ,0) = C Ly (2.21) 

n Zl tor(Ly, 1) = (^J ~lj (2-22) 

n z ,tor( L y, d) = [J^^ rf J for 2 < d < L y (2.23) 
where C n is the Catalan number which occurs in combinatorics and is defined by 

C n = -±-( 2n \. (2.24) 



(n + 1) \ n 

The first few Catalan numbers are C\ = 1, G% — 2, C3 — 5, C4 = 14, and C5 = 42. Special 
cases of n Z) tor(L yi d) that are of interest here include 

n z ,tor(L y ,L y ) = 1 (2.25) 

nz,tor{L y , Ly - 1) = 2L y for L y > 3 (2.26) 

In Table H] we list the first few numbers nz,tar{L y ,d). We notice the following recursion 
relations: 

nz,tor(L y , 1) = n z ,tor(L y , 0) - n z ,tar(L y - 1, 0) + 2n z ,tor(L y -1,1) + n z ,tor(L y - 1, 2) 
for L y > 2 

nz,tor{L y , 2) = -n z ,tor{L y - 1, 0) 4- 2n z ,tor{L y - 1, 1) + 2n z ,tor{L y - 1, 2) 4- n z ,tor{L y - 1, 3) 
for L y > 2 
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nz,tor(L y , d) = n z ,tor(L y -l,d-l) + 2n z<tor (L y -l,d) + n Z)tor (L y - 1, d + 1) 

for L y > 3, 3 < d < L y . (2.27) 

The last line here is the same as the recursion relation for n Z)Cyc (L y , d) when 1 < d < L y for 
cyclic strips (which was denoted as n z (L y ,d) in By substituting q = into eq. (j2.19|) 

and using eq. (j2.15|) . we have the relation 

Ly 

^2n z ,tar{L y ,d)(-l) d =0 for A = sq, tri, he . (2.28) 

d=0 



TABLE I: Table of numbers riz,tor(L y ,d). See text for general formulas. 
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6 


1 
















4 


14 


35 


28 


8 


1 














5 


42 


126 


120 


45 


10 


1 












6 


132 


462 


495 


220 


66 


12 


1 










7 


429 


1716 


2002 


1001 


364 


91 


14 


1 








8 


1430 


6435 


8008 


4368 


1820 


560 


120 


16 


1 






9 


4862 


24310 


31824 


18564 


8568 


3060 


816 


153 


18 


1 




10 


16796 


92378 


125970 


77520 


38760 


15504 


4845 


1140 


190 


20 


1 



Another way to realize the number n ztor (L y ,d) is as follows. For d = 0, i.e. no color 
assignment, the partitions are identical to those for cyclic strips, and the number of these 
partitions is n Zit ar(L y ,0) = n Zcyc (L y ,0) = C^ y given by eq. (|2.21|) . Recall that the number 
of partitions of size n with k components is g iven by the Narayana number N(n, k), which 
is sequence A001263 in [3^ (see also Ref. 37] and the references therein), given by 
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The first few rows of the triangle of Narayana numbers, or Catalan triangle, are shown in Fig. 
0] It is clear that the sum of terms in each row in the triangle is equal to the corresponding 
Catalan number: 




= n z<tor (L y ,0) . (2.30) 



For d = 1, we need to assign a color for each component of the level partitions. Therefore, 
we have 

n zM ) - g k N { L„ *) = £ft:}) (/_\) = gCi) 

1 {2L V \ L v + 1 L v + 1 . . . 

= 2 v l j = ~V" = - Y - "* 1 "-^' 0) • (2 - 31) 

Let us compare nz,tar{L y ,d) with the corresponding number nz, cy c(L y ,d) given in Table 3 
of In general, n ztor {L y ,d) is larger than n z , C yc(L y ,d), except for d = where they 

are the same. This is a consequence of the fact that certain forbidden partitions for cyclic 
strips for d > 1 are allowed for toroidal strips with rotational translation. We list the first 
few nz,tor(L y , d) — nz, C yc{L y , d) in Table ITT1 Comparing Table HP with Table HI we infer the 
relation 

n z ,tor{L y , d) - n z , C yc{L y , d) = n z>tor (L y , d + 1) for d > 2 . (2.32) 

This can be understood as follows. The extra partitions for toroidal strips, as contrast 
with cyclic strips, are those with removable crossings. That is, for any connection between 
two non- adjacent white circles, all the color assignments should be inside or outside the 
connection while these two white circles can be either color- assigned or not. Therefore, there 
are three possible correspondences from the partitions for n Zjt0 r(L y , d + 1) to the partitions 
for nz,tor(Ly,d) — n Zcyc (L y ,d) with d > 2: (i) if the top color-assigned vertex and the 
bottom color-assigned vertex are not within a connection, then connect these two vertices 
(so that only one color is assigned to them); (ii) if all color-assigned vertices are within a 
connection and the connection is also assigned a color, then remove the assignment for the 
connection; (iii) if all color-assigned vertices are within a connection and the connection is 
not assigned a color, then connect the top and bottom color- assigned vertices (within that 
connection). These account for all of the extra partitions that toroidal strips have, relative 
to cyclic strips. As an example, consider the case of width L y = 4 and degree d = 2 in eq. 
(J2.18|) . For this case there are eight partitions that are allowed for toroidal strips but not 
for cyclic strips. These are 13,2; 14,3; 14,2; 14,2,3; 24,3; 14,23; 124,3; 134,2. Here we 
only use over line to indicate the color assignment, since permutations of black circles are 
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not considered for nz,tor(L y , d). These partitions are in one-to-one correspondence with the 
partitions for L y = 4 and d = 3 (without color permutation), namely 1,2,3; 1,3,4; 1,2,4; 
14, 2, 3; 2, 3, 4; 23, 1, 4; 12, 3, 4; 34, 1, 2. Eq. (I2~!a is equivalent to the relation 

Ly 

n z ,tor(L y , d) = } j n ZjCyc (L y , d') for d > 2 . (2.33) 
d'=<f 

1 

1 1 
1 3 1 

16 6 1 

1 10 20 10 1 

1 15 50 50 15 1 

1 21 105 175 105 21 1 
1 28 196 490 490 196 28 1 
1 36 336 1176 1764 1176 336 36 1 

FIG. 4: Triangle of Narayana numbers. 

The construction of the transfer matrix for each level (= degree) d can be carried out 
by methods similar to those for cyclic strips 23]. Notice that for the honeycomb lattice, 
the number of vertices in the transverse direction, L y , should be an even number, and the 
smallest value without degeneracy is L y = 4. Using the bases described above (e.g. eqs. 
()2.16|) - (j2.18J) for 2 < L y < 4), we define JL y ,d,i,i+\ as the join operator between vertices i and 
i + 1 and Dh y 4,i as the detach operator on vertex i for each subspace d. As compared with 
the situation for cyclic strips, the corresponding toroidal strips with the same L y and d have 
one more join operator, namely J Ly ,d,L y ,i = JL y ,d,L y ,L y +i- The transfer matrix T ZAjLy4 for 
each d is the product of the transverse and longitudinal parts, Hz,\,L v ,d and Vz,\,L y ,di which 
can be expressed as 

Ly 

i=l 

Ly 

Hz,tri,L y ,d = J L y +\,d,Ly+l,l + v JLy+l,d,i,i+l) 

1=1 
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TABLE II: Table of numbers nz,tor(L y , d) — nz tC yc(L y , d). 



J-* y \a 


n 

u 


1 

_L 


9 




4 


K 


U 


7 


a 

O 


q 


1 


1 


u 


u 




















o 

z 


U 


U 


U 


















Q 

o 


U 


1 


1 



















4 


U 


I 


o 


1 

















5 


o 


36 


45 


10 


1 


o 












6 





165 


220 


66 


12 


1 













7 





715 


1001 


364 


91 


14 


1 











8 





3003 


4368 


1820 


560 


120 


16 


1 









9 





12376 


18564 


8568 


3060 


816 


153 


18 


1 







10 





50388 


77520 


38760 


15504 


4845 


1140 


190 


20 


1 






Ly/2 Ly/2 

Hz,hc,L y ,d,l = Y\_ (I v ^Ly,d,2i-l,2i) , Hz,hc,L y ,d,2 = JJ (/ + V J L y ,d,2i,2i+l) 

i=l i=l 

Ly 

Vz,sq,L y ,d = Vz,hc,L y ,d = J^C 1 ^ + ^i s ,d,i) 

i=l 

Ly 

^iri.L^d = ^+1,^,1 [][(/ + ^Lj+MAi+i)^ + -D^+iAi+i)] » (2-34) 

where [z/] denotes the integral part of u, and / is the identity matrix. Notice that for 
triangular strips with periodic transverse boundary conditions, one has to work with width- 
(L y + 1) matrices and then identify the two end vertices, as shown in 0|. We have 

Tz,sq,L y ,d ^Z,sq,L y ,dH Z,sq,L y ,d > Tz t tri,L y ,d ^Z,tri,L y ,dH Z,tri,L y ,d 

Tz,hc,L y ,d = (Vz,hc,L y ,dH Z,hc,Ly,d,2)(Vz,hc,L y ,dH z,hc,L y ,d,l) = Tz,hc,L y ,d,2Tz,hc,LyA^ ' (2-35) 

As discussed above, the definition of does not include the possible permutations of 
black circles. In principle, there can be d\ eigenvalues with their coefficients bj added up to 
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. In our explicit calculation for toroidal strips, decomposes into two b^ (with possible 
multiplication by an even integer) for the square and honeycomb lattices when 1 < d < L y 
up to L y = 4 and for the triangular lattice when 1 < d < L y up to L y = 3. As discussed 
above, we know that there is only one coefficient for both d = and d — 1, i.e., b^ = 1, 
ftC 1 ) — q — 1. For d = 2, 3, the coefficients are 



b? = 




-3) 




u 2 — 






-2) 


bf = 


|(«- 


1)(9 2 


- 5g + 3) 










6 (3) _ 


1 , 
6 9(9 


-2)(? 


-4) 



(2.36) 

with possible multiplication. In contrast to cyclic and Mobius strips, certain eigenvalues for 
toroidal strips may appear in more than one level, so that their coefficients are the summation 
of the corresponding fe^'s. 

For the Klein bottle strip of the square and triangular lattices or the honeycomb lattice 
with L y even, the sum of coefficients is the same as the sum for the corresponding strip with 
Mobius longitudinal boundary conditions, 

_ "W" f <f,l* for even L, 

CZ ' L - Kb= U CZ - L - KbJ= \ q i^ for odd L, ' <2 ' 37> 

Consider the coefficients for Klein bottle strips. For the square and honeycomb lattices, when 
the longitudinal boundary condition is changed from toroidal to Klein bottle, we observe the 
following changes of coefficients: 

6 (o) ^ ±6 (o) 



6 (i) ^ ±6 (i) 

&S 2) -> ±b? 

bf - ±bf 
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(2.38) 



We find that certain coefficients become zero when the boundary condition is changed from 
toroidal to Klein bottle, and therefore the number of eigenvalues for Klein bottle strips 
always appears to be less than the number for the corresponding toroidal strips. This has 
been observed before in Ref . 0, |^ . The partition functions for Klein bottle strips have 
the same expressions as eqs. (|1.5|) and (jl,15|) with replaced by appropriate coefficients, 
as in eq. (|2.38j) . Illustrative calculations will be given below. 



III. ZERO-TEMPERATURE POTTS ANTIFERROMAGNET 

For the chromatic polynomials, it is necessary that adjacent vertices are not assigned the 
same color, i.e., in the partition diagrams, adjacent vertices cannot be connected by color- 
connection lines. For the square and triangular lattices, a transverse slice is a circuit graph; 
while for the honeycomb lattice, a transverse slice is a set of L y /2 two- vertex tree graphs. 



A. Square and triangular lattices 

For the square and triangular lattices, the second partition of the degree d = subspace 
and the third partition of the degree d = 1 subspace for width L y = 2 in Fig. El are not 
allowed. For L y = 3 in Fig. we keep the first partition in the d = subspace, the first 
three partitions of the d = 1 subspace, and the first six partitions of the d = 2 subspace. 
Let us denote the reduced size of the transfer matrix for the square and triangular lattices 
as np t tor(L y ,d), excluding permutations of black circles. The sum of all coefficients is the 
dimension of the total transfer matrix (independent of L y ), which is the chromatic polynomial 
of the circuit graph given by Theorem 10 of [3^, so that for L y > 2, 

Ly 

Cp,A, Ly = J2 n P,tor(L y ,d)b {d) = (q- l) Ly + {q- l)(-l) L y for A = sq,tri . (3.1) 

d=0 

By substituting the expression for £>w in eq. (|2.19|) . we determine the n Pjtor (L y , d). We 
list the first few numbers n Pjtor (L y ,d) in Table ITTTl The strip with L y = 1 is obtained by 
identifying the free boundaries of the cyclic strip with L y = 2, namely, each vertex is attached 
by a loop such that its chromatic polynomial is identically zero. The column npj or {L y , 0) in 
Table ITTTl is the number of non-crossing non-nearest-neighbor partitions of L y vertices on a 
circle, and is given by the Riordan number Rl |4l|. If we compare np ttor (L y , 0) with the 
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correspondin g nu mber np )Cyc (L y ,0) (which was denoted as np(L y ,d)) for cyclic strips, given 
in Table 1 of Li5L we find the relation 



np, C yc{Ly, 0) - n P j or (L y , 0) = rip !tor (L y - 1, 0) for L y > 2 . (3.2) 

This can be understood as follows. For toroidal strips we have additional bonds connecting 
the top and bottom vertices of the corresponding cyclic strips, so these vertices cannot have 
the same color. That is, the d = partitions for cyclic strips with the top and bottom 
vertices connected (having the same color) is not allowed for the corresponding toroidal 
strips. The number of these partitions is given by n Ptor (L y — 1,0). This is clear, since 
we know np iCyc (L y ,0) = Ml \ |35j|, where M n is the Motzkin number, and the relation 
M n = R n + R n+ \. We notice the following recursion relations: 

n P)tor (L y , 1) = n P j or (L y , 0) + n P j or (L y -1,1)+ n P)tor (L y - 1, 2) + (-l) Ly for L y > 2 
n Pjtor (L y , 2) = -n Pitor (L y - 1, 0) + 2n Pttor (L y -1,1) + n P)tor (L y -1,2) 

+n P)tor (L y -1,3) + (-l) L y for L y > 2 
np,tor(L yj d) = n P j or (L y - 1, d - 1) + n Pitor (L y -l,d) + n P j or (L y — l,d+l) 

for L y > 3, 3<d<L y . (3.3) 

The last line here is the same as the recursion relation for np )Cyc (L yi d) when 1 < d < L y for 
cyclic strips [35]. By substituting q = into eq. (|3.1|) and using eq. ()2.15|) . we have 

Ly 

^2n P<tor (L y ,d)(-l) d = forA = sg,tri. (3.4) 

d=0 

Let us denote the number of diagonal dissections of a convex n-gon into k regions as 
D(n, k), which is sequence A033282 in (3^. (See also Ref. |3t]].) 

D(n,k) = Tf?~ 3 V n | fc ~ 2 l iorn>3,l<k<n-2. (3.5) 



k\k — lj\ k - 1 

The first few numbers are £>(3, 1) = 1, D(4, 1) = 1, D(4, 2) = 2, D(5, 1) = 1, D(5, 2) = 5, 
D(5, 3) = 5. We also know that the Riordan number for n > 1 is given by the number of 
dissections of a convex polygon by nonintersecting diagonals with n + 1 edges. To relate 
D(n, k) with the Riordan number, define D'(n, k) as 

D'(n,k) = D(n + 2- k,k) = \( U ~ 1 ~ ) ( , " ) for n > 2 . (3.6) 

k\ k — 1 \k — 1 



19 



TABLE III: Table of numbers np t tor{L y ,d) for the square and triangular lattices. 



L \ H 


n 


i 
± 


9 


o 


A 
-± 


x. 

O 


U 


7 


8 
o 


Q 


1 


-1 

1 





U 




















o 
z 


1 


o 
Z 


1 


















3 


1 


3 


3 


1 
















A 

4 


3 


1U 


10 


4 


1 














5 


6 


25 


30 


15 


5 


1 












6 


15 


71 


90 


50 


21 


6 


1 










7 


36 


196 


266 


161 


77 


28 


7 


1 








8 


91 


554 


784 


504 


266 


112 


36 


8 


1 






9 


232 


1569 


2304 


1554 


882 


414 


156 


45 


9 


1 




10 


603 


4477 


6765 


4740 


2850 


1452 


615 


210 


55 


10 


1 



It is clear that the entries of D'(n,k) are zero for [n/2] + 1 < k < n. Let us also define 
D'(l, 1) = = rip itor (l, 0), and list the first few rows of the triangle of D'(n, k) for n > 1, 1 < 
A; < n in Fig. The summation of each row in the triangle gives the Riordan number: 

[Ly/2] 

D'(L y , k) = R Ly = n P>tor {L y , 0) . (3.7) 

k=l 

For d = 1, we need to assign a color for each component of level partitions, so that 

[Ly/2] 

n Pt tor(L y , 1) = (Ly + 1 - k)D'(L y , k) = 

k=l 

For entries np^cyr{L y ,d) in Table ILTTl with d > 1, we find the following relation with 
known number sequences. Motivated by eq. ()3.2|) . let us define n' Ptor (L y ,d) such that 
n Ptor (l,d) = n P)tor {l,d) and 

n p,tor(A/> rf ) = n P,tor(L y , d) + n PM .(L y - 1, d) for L y > 2 . (3.9) 

We list the first few numbers n' ptor .(L^d) in Table ITVl Comparing these with the triangular 
array given by sequence A025564 in |36( (which result from pairwise sums of entries in the 




20 





1 
1 
1 2 
1 5 
1 9 5 
1 14 21 
1 20 56 14 
1 27 120 84 

FIG. 5: Triangle of D'(n, k). 



trinomial array (sequence A027907)), we find that n' Ptor (L y ,d) for d > 2 are columns in 
the array. The first few rows of this triangular array are listed in Fig. |H1 Let us denote a 
particular entry as T(n, fc); then T(n, k) = T(n — 1, k — 2) + T(n — 1, k — 1) + T(n — 1, k) 
starting with rows [1],[1, 2, 1],[1, 3, 4, 3, 1]. We find the relation 



n 



Ptor {L y , d) = T(L y + l,L y + l - d) = T(L y + 1, Ly + 1 + d) forrf>2. (3.10) 



We list the first few n' ptor (L y ,d) — np jCyc (L y ,d) for L y > 2 in Table M Because of eq. 
we have n' ptor (L y , 0) - n PiCyc (L^, 0) = 0. For d = 1, we find n' ptor (L y , 1) - n PjCyc (L y , 1) 
is given by sequence A005775 in [36], the elements of which are the numbers of compact 
rooted directed animals of size L y > 3 having 3 source points. Combining the results in 
Table Eland Table HVl we infer the relation 



n P,tor( L w d ) ~ n P,c V c{L y , d) = n' P>tor (L y , d+1) for d > 2 , (3.11) 

or equivalently, 

Ly 

n 'p,tor( L yi d) = n PtCyc (L y , d') for d > 2 . (3.12) 

We recall that the sum of the coefficients in (|1.9|) for a Klein bottle strip of the square or 
triangular lattice of width L y is (independent of L x ), given as Theorem 11 in j^], is 

Cp,L y ,Kb = Cp,L y ,Kb,j = (3.13) 
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TABLE IV: Table of numbers n' Ptor (L y ,d). 



L y \ a 


u 


i 
i 


o 
Z 


q 


A 

4 





D 


7 


o 
o 


y 


1 n 

1U 


1 


























2 


1 


2 


1 


















3 


2 


5 


4 


1 
















4 


4 


13 


13 


5 


1 

















y 






1 Q 

iy 


ft 



1 
1 












6 


21 


96 


120 


65 


26 


7 


1 










7 


51 


267 


356 


211 


98 


34 


8 


1 








8 


127 


750 


1050 


665 


343 


140 


43 


9 


1 






9 


323 


2123 


3088 


2058 


1148 


526 


192 


53 


10 


1 




10 


835 


6046 


9069 


6294 


3732 


1866 


771 


255 


64 


11 


1 



B. Honeycomb lattice 

For a strip of the honeycomb lattice with toroidal boundary conditions, the width L y 
must be even. In each transverse slice of such a strip, the L y vertices are connected in a 
pairwise manner, just as for the strips with cyclic boundary conditions. Let us denote the 
reduced size of the transfer matrix as n Pjtorhc (L y ,d), without considering permutations of 
black circles in the partition diagram. The sum of all coefficients is the dimension of the 
total transfer matrix (independent of L x ), which is the same as eq. (6.19) of [3^ for cyclic 
strips, namely 

Ly 

C P , hc , Ly = J2 n P,tor,hc(L y ,d)b^ = {q{q - 1)) L * /2 . (3.14) 

d=0 

By substituting the expression for in eq. (|2.19J) . we determine the np t tor,hc(Ly, d). We list 
the first few numbers np )tor ^ c {L yi d) in Table IV II We infer the following recursion relations 
for L y > 4: 

np,tor,hc(Lyi 0) = —np ttar ,he{Ly — 2,0)+ ^np tor hc {L y — 2, 1) — np^ r.hSJ J y ~ 2, 2) 

—np t tor,hc(.Ly ~ 2, 3) 
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1 

1 3 
14 8 
1 5 13 22 
1 6 19 40 61 
1 7 26 65 120 171 
1 8 34 98 211 356 483 
1 9 43 140 343 665 1050 1373 



1 

2 1 

4 3 1 

10 8 4 1 

26 22 13 5 1 

70 61 40 19 6 1 

192 171 120 65 26 7 1 

534 483 356 211 98 34 8 1 

1500 1373 1050 665 343 140 43 9 1 



FIG. 6: Triangle of T{n,k). 



TABLE V: Table of numbers n' Ptor (L y , d) — np^ cyc (L y7 d). 



L y \ d 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


2 



























3 





1 


1 



















4 





4 


5 


1 

















5 





14 


19 


6 


1 















6 





45 


65 


26 


7 


1 













7 





140 


211 


98 


34 


8 


1 











8 





427 


665 


343 


140 


43 


9 


1 









9 





1288 


2058 


1148 


526 


192 


53 


10 


1 







10 





3858 


6294 


3732 


1866 


771 


255 


64 


11 


1 






np,tor,hc(Ly, 1) — 
np,tor t hc{L y , 2) = 



— 3np jtori / lc (Lj / — 2, 0) + 10np ttorthc (L y — 2, 1) + np )t0 r,hc(L y — 2, 2) 
— 3np itori / lc (L ?/ — 2, 0) + Sn Pttorjhc (L y — 2,1)+ ^n P ^ tor . thc (L y — 2, 2) 
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J r^ n P ! tor,hc(-^'y — 2, 3) + np t tar,hc\Ly ~ 2, 4) 

° n 'P,tor,hc^J-'yi 3) = ~ n P,tor,hc(L y — 2,0) + 2np tor } lc (L y — 2, 1) + Sripjor ^Ly — 2, 2) 

+^ n P,tor,hc(L y — 2, 3) + 3np jtor!hc (Ly — 2,4)+ n P tor hc (L y — 2, 5) 

np,tor,hc(-L y , d) = np jt0 r,hc(L y — 2, d — 2) + ^np t tor,hc(L y — 2, d — 1) + 4rip )tor) /i C (L 2/ — 2, d) 

+3np t tor,hc{Ly -2,d + l) + n PM ;hc( L y - 2, d + 2) for 4 < d < L y (3.15) 
By substituting g = into eq. (|3.14j) and using eq. ()2.15|) . we have 

£)np, tor , hc (L 1 „tO(-l) <i = 0. (3.16) 

d=0 



TABLE VI: Tabic of numbers np t tor,hc(L v ,d) for the honeycomb lattice. 



L y \ d 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


2 


1 


2 


1 


















4 


6 


18 


17 


6 


1 














6 


43 


179 


213 


108 


39 


9 


1 










8 


352 


1874 


2518 


1512 


721 


264 


70 


12 


1 






10 


3114 


20202 


29265 


19425 


10800 


4953 


1830 


525 


110 


15 


1 



It is clear that the number of degree d = partitions np 7 tor,hc(-Ly>Q) is the same as the 
corresponding number n PtCyC)hc (L y , 0) for the corresponding cyclic strips (which was denoted 
as np(hc, L y ,d) in 0]). Combining these results with those in Table 2 of 0], we observe 
that 

It i \ P,cyc,hc,L y ,X , n 1 _s 

np,tor,hc( L y> !) = 2 ^ ^ 

where Np tCyCt h c ,L y ,x is the total number of A for cyclic strips. We list the first few 
np t tor,hc(.L y ,d) — np tCyCt hc(Ly,d) in Table IVIll Compare Table IVlTl with Table ED it is easy 
to see the relation 

np, tor ,hc( L y> d) - np iC y Cthc (L y , d) = n PttorM (L y , d + 1) for d > 2 , (3.18) 
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or equivalently, 
We also have 

T r l'P t tor,hc\^ J y> 1) 



np,tor,hc( L yi d) = ^2 n P,cyc,hc(L y , df) for d > 2 . 



(3.19) 



d'=d 



np t cyc,hc(Ly,l) — ^^( — 1)^ +ln P,tor,hc(L y , d!) — 71 P !C y C ,hc(L y , d!) . 

(3.20) 



d'=3 



odd d'=3 



TABLE VII: Table of numbers npjor,hc(L y ,d) — np iCVCi hc{L y , d) for the honeycomb lattice. 



L y \ d 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


2 



























4 





5 


6 


1 

















6 





77 


108 


39 


9 


1 













8 





996 


1512 
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264 
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1 
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12177 


19425 


10800 


4953 


1830 


525 


110 


15 


1 






IV. PROPERTIES OF TRANSFER MATRICES AT SPECIAL VALUES OF 

PARAMETERS 

In this section we derive some properties of the transfer matrices Tz,\,L y ,d at special values 
of q and v, and, correspondingly, properties of TT,\,L y ,d at special values of x and y. 



A. v = 

From (jl.lj) or (jl.lUJI it follows that for any graph G, the Potts model partition function 
Z(G, q, v) satisfies 

Z(G,q,0) = q n ^ . (4.1) 

Since this holds for arbitrary values of q, in the context of the lattice strips considered here, 
it implies that 

(T z ,A,L y> d)v=o = for 1 < d < L y , (4.2) 
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i.e. these are zero matrices. Secondly, restricting to toroidal strips for simplicity, and using 

the basic results n = L y L x = L y m for A = sq, tri and n = 2L y m for A = he, eq. (|4.1|) 
implies that 

{ L y m f or _ g _ ^ 
q IL v m for A = he 

With our explicit calculations, we have 

{Tz,!v,L y ,o) jk = for v = 0, and j > 2 , (4.4) 

i.e., all rows of these matrices except the first vanish, and 

{T z ,A,L y ,a)u = q pLy for v = 0, (4.5) 

where p was given in eq. ()1.18|) . For this v = case, since all of the rows except the first 
are zero, the elements (Tz,A,L y ,o)ik for > 2 do not enter into Tr[(Tz,A,L v ,o) m ), which just 
reduces to the m'th power of the (1,1) element: 

Tr[(T ZA>Lyfi r] = [(T z>A ,L v ,o)iir ■ (4.6) 

Corresponding to this, all of the eigenvalues Xz,A,L y ,d,j vanish except for one, which is equal 
to (7^^,^,0)11 in eq. (|4.5j) . As will be seen, this is reflected in the property that det(T Z) A,L y ,d) 
has a nonzero power of v as a factor for all of the lattice-A strips considered here. Note that 
the condition v — is equivalent to the Tutte variable condition y — 1. 



B. q = 

Another fundamental relation that follows, e.g., by setting q = in eq. (jl.lOj) . is 

Z(G,0,v) = 0. (4.7) 

Now the coefficients &w evaluated at q = satisfy eq. ([2.15|) . Although b^ d > decomposes into 
two bj for d > 1, one of them is equal to zero and the other is (— l) d at q = for the cases 
we have, namely = 0, 63 = 1, 6^ = — 1, 62 = 0- Hence, in terms of transfer matrices, 
we derive the sum rule 

E E( A w,,<*,iM) m - E E( A w»^(o,v)) m = . (4.8) 

CKeven d<L y j l<odd d<L y j 

This is similar to a sum rule that we obtained in Ref. j^. Since it applies for arbitrary m, 
it implies that there must be a pairwise cancellation between various eigenvalues in different 
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degree-<i subspaces, which, in turn, implies that at q = there are equalities between these 
eigenvalues. For example, for the L y = 2 strip of the square lattice and q = 0, two of the 
eigenvalues of Tz )a q,2,i become equal to the eigenvalues of Tz, aqj 2,o, while the third eigenvalue 
of Tz, sq ,2,i becomes equal to \z,aq,2,2 = v 2 . Note that setting q = does not, in general, 
lead to any additional vanishing eigenvalues for the Tz,K,L y ,d of the A = sq, he lattices, and 
hence our formulas below for Tz,A,,L y ,d do not contain overall factors of q. However, for the 
triangular lattice, certain eigenvalues do vanish at q — 0. 

In terms of the variables x and y in the Tutte polynomial, the value q = is equivalent 
to the value x = 1 (unless v = 0). In contrast to the vanishing of Z(G,q,v) at q = 0, the 
Tutte polynomial T(G, l,y) is nonzero for general y. This different behavior can be traced 
to the feature that in eq. 1)1. 14j) . Z(G,q,v) is proportional to T(G,x,y) multiplied by the 
factor (x - l) k ^ = (x - 1), so at x = 1, Z(G, 0, v) = even if T(G, 1, y) ^ 0. 

C. v = —1, i.e., y = 

As discussed above, the special value v = —1, i.e., y — 0, corresponds to the zero- 
temperature Potts antiferromagnet, and in this case the Potts model partition function 
reduces to the chromatic polynomial, as indicated in eq. (jl.6|) . In the previous two sec- 
tions we have determined how the dimensions riz t tor{L y ,d) for the square, triangular, and 
honeycomb lattice strips reduce to the dimensions n P)tor (L yi d). In all cases, in each degree- 
d subspace for < d < L y — 1 for A = sq,tri,hc, some eigenvalues vanish, so that 
np t tor,A(L y ,d) < nz,tor,x{Ly,d) and hence det(Tz t x,L y ,d) — at v = —1. This property is 
reflected in the powers of (v + 1) and y that appear, respectively, in our formulas below for 
det(T z ,A,L y ,d) and det(T T ,A,L y ,d)- 

D. v = —q, i.e., x = 

For the graph G = G(V, E), setting x = and y = 1 — q in the Tutte polynomial T(G, x, y) 
yields the flow polynomial F(G,q), which counts the number of nowhere-0 g-flows (without 
sinks or sources) that there are on G [43]: 

F(G, q) = (-l) e(G) ~ n(G)+1 T(G, 0, 1 - q) . (4.9) 

Therefore, the flow polynomial for toroidal lattice strips has the form 

Ly 

F(A,L y x L x ,tor.,q) = J2 b f ^AMiT ■ ( 4 -!0) 

d=0 j 



27 



We found that for all L y > 1 and < d < L y — 1 for A = sq, he, when one sets x = in the 
Tutte polynomial, or equivalently, q = —v in the Potts model partition function, some of 
the eigenvalues in each degree-d subspace vanish, and hence det{Tz,k,L y ,d) and det{TT^L y ,d) 
vanish. This is reflected in the powers of (q + v) and x that appear, respectively, in our 
formulas below for det{Tz,k,L y ,d) and det{TT^L y ,d)- The expressions for the determinants for 
strips of the triangular lattice are more complicated because of the fact alluded to above 
that one has to work with width- (L y + 1) matrices and then identify the two end vertices. 



V. GENERAL RESULTS FOR TOROIDAL STRIPS 

In this section we present general results that we have obtained for transfer matrices of 
toroidal strips and their properties. These are valid for arbitrarily large strip widths L y > 2 
(as well as arbitrarily great lengths). 



A. Determinants 

For the square lattice, we find 

det(T T;Sq>Ly4 ) = (xyf W^-M) ( 51 ) 

where n z ,tor(L y , d) was given by eq. (j2.21J) - ()2.23|) . This result applies for all d, i.e., < d < L y 
since nz,tor(L y — 1, d) = for d > L y — 1 (c.f. eq. ()2.20|) ). This is equivalent to the somewhat 
more complicated expression for det(Tz )S q,L y ,d)' 



det(T z ^ Ly4 ) = v dl L y^oAL y> d) M + — j (1 + t>) . (5.2 



Q 



v 



d\ L y n z ,t or (Ly-l,d) 



These determinant formulas can be explained as follows. By the arrangement of the par- 
titions, as shown in Figs. |2]and|21 the matrix JL y ,d,i,i+i has the lower triangular form and 
the matrix -DL y ,d,i has the upper triangular form. From the definition of the transfer ma- 
trix in eq. (j2.r>5j) . the determinant of Tz, S q,L y ,d is the product of the diagonal elements of 
H z ,sq,L v ,d and V z , S q,L v ,d- The diagonal elements of H z , S q,L y ,d have the form (1 + v ) r , where r 
is the number of edges in the corresponding partition, and the diagonal elements of Vz, S q,L y ,d 
have the form v Ly (l + q/v) s , where s is the number of vertices which do not connect to any 
other vertex in the corresponding partition. Therefore, the power of (1 + v) in eq. ()5.2j) is 
the sum of the number of nearest-neighbor edges of all the (L y , d) -partitions. Let us com- 
pare the (L y , d)-partitions and [L y — l,d) partitions. Since each edge of a (L y , d) -partition 
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corresponds to adding a new vertex in a (L y — 1, d)-partition (at L y possible places) and 
connecting it with the vertex above it, the power of (1 + v) is d\ L y n ztor (L y — 1, d), where 
the factor d\ comes from the permutation of black circles. It is clear that the power of 
v is d\ L y nz,tor(L y ,d). The power of (1 + q/v) in eq. ()5.2|) is the sum of the number of 
unconnected vertices of all the (L y , (i)-partitions. Now consider the (L y , <i)-partitions and 
(L y — l,d) partitions again. Since each unconnected vertex of a (L y , d)-partition corresponds 
to adding a unconnected vertex to a [L y — 1, d)-partition in L y possible ways, the power of 
(1 + q/v) is again d\ L y n Zttor (L y - 1, d). 

Next, taking into account that the generalized multiplicity is U- d \ we have, for the total 
determinant 

Ly 

L v n z , tor (L y ,d)b<- d 1 1 mi \Lyn z , tor (Ly-l,d)b^ 



det{T z ^ Ly ) = Y[{y-l) L * nz ^ L ^ bid \xy) L * 



d=0 



Using eq. together with eq. (l2~271|) so that Y^dU n ZMr (L y -l, rf)6 (d) =Eto 1 n zM L y- 

l,d)b( d \ we have, finally, 



det(T z>sqjLy ) = (y - l) L ^ v {xyf^ 1 . (5.4) 



This agrees with the conjecture given as eq. (3.70) of our earlier Ref. |8( for the determinant 
of the transfer matrix of the toroidal strip of the square lattice with arbitrary width L y . 
For the honeycomb lattice, we find 

det(T Wy4 ) = (x 2 yf W^-M) . (5 5) 

Equivalent ly, 

r/ n \ 2 -i d\ L y nz tor(L y — l,d) 

det(T ZM , Ly4 ) = (v 2 f Wo,(^.rf) + £j (i + v )j . (5. 6 ) 

This can be understood as follows: by an argument similar to that given before, the power 
of (1 + v) is the same as for the square lattice case. Comparing Tz )S q,L y ,d an d T z ^c,L y ,d in 
eq. (|2.H5j) . one sees that V z ^ c ,L y ,d = Vz, S q,L y ,d has been multiplied twice for the honeycomb 
lattice, so that the powers of v and (1 + q/v) become twice of the corresponding powers for 
the square lattice. 

Taking into account that the generalized multiplicity is b^ d \ the total determinant for the 
he lattice is given by 

det(T T>hc r) = (x'y)^- 1 . (5.7) 
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Ly-1 



Ly-1 



1 % 

2 



Ly-1 



J L v 



FIG. 7: m = 1 graph for the toroidal square lattice. 



Equivalent ly, 



det(T : 



,hc.L v 



,2\L y q L v 



1 + 



(5.f 



It is clear that det(T TjhCjLyjd ) is related to det(T TsqLyjd ) by the replacement x — > x 2 (hold- 
ing y fixed). This, together with the fact that nz,tor{L y , d) is the same for these lattices means 
that the total determinants det(T T>h ^ Ly ) is related to det(T T;Sq>L ) by the same respective re- 
placements. Correspondingly, det(Tz,hc,L y ,d) is related to det{Tz, S q,L y ,d) by the replacements 
of the respective factors v by v 2 (cf. eq. (11.18(0 and (1 + q/v) by (1 + q/v) 2 . 



B. Traces 

The trace of the total transfer matrix is the m = 1 case in eq. (jl.5j) or ()1.15|) . For strips of 
the square lattice, this corresponds to a L^-vertex circuit with a loop attached to each vertex, 
as illustrated in Fig. [7| We have Tr{TT yS q,L y ) given by the corresponding Tutte polynomial, 

Tr{T My )=y L y(^-l + ^^j . (5.9) 

Equivalent ly, we find 

Tr(T z ^ Ly ) = (1 + v) L *[(v + qfy + (q - l)v L y] . (5.10) 

In principle, we can also consider the m = 1 case for the Klein bottle strips, but the result 
is not as simple as that listed here. 

For the triangular lattice, the total trace corresponds to a L y -vertex circuit with each 
edge doubled and with a loop attached to each vertex as illustrated in Fig. |H1 Therefore, 
the trace is given by the Tutte polynomial of this graph 

Tr(T T>triiLy ) = y L y \^y - l)(y + 1) L - + £j (y + iy(x + . (5.11) 
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Ly - 1 Ly-1 




FIG. 8: m = 1 graph for the toroidal triangular lattice. 



Ly-1 



1 <fZ3^ 
2 



4 C3 



Ly — 1 Ly — 1 < t 

Ly Ly j tl S 



FIG. 9: m = 1 graph for the toroidal honeycomb lattice with even L y . 



Equivalent ly, 



Tr(T ZM . Ly ) = q(v + 1)^ ( {v 2 + 2v) L * + ^ {v 2 + 2 V y{v 2 + 2v + g) 1 *" 1 

j=0 



(5.12) 



For the honeycomb lattice, the total trace corresponds to a 2L y -vertex circuit with every 
other edge doubled as shown in Fig. El Hence, the trace is given by the Tutte polynomial of 
this graph 

^^-j +(y-l)(y + l) L »+J2(y + iy(x + y) L »- 1 -i . (5.13) 

Equivalently, 

Tr(T ZM;Lv ) = {v 2 + 2v + q) L y ({v + q) L y -v L y) 



+qv Ly (v 2 



+ 2^) L " + (v 2 + 2w ) V + 2w + g^ -1 "' ] • (5-14) 

3=0 
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C. Eigenvalues for d = L y and d = L y — 1 for A = sq, he 

It was shown earlier that the A's for a strip with Klein bottle boundary conditions 

are a subset of the A's for the same strip with torus boundary conditions. From eq. (|2.23|) 
one knows that there is only one A, denoted as \z,A,L y ,L y , for degree d = L y for strips of the 
square and honeycomb lattices because all the permutation of black circles are equivalent. 
That is, for this value of d, the transfer matrix reduces to 1 x 1, i.e. is a scalar. We found 
that for a toroidal or Klein bottle strip of the square, or honeycomb lattice with width L y , 

^T,A,L y ,L y = 1 ■ (5.15) 
Equivalently, in terms of Potts model variables, 

\z,s q ,L y ,L y = v Ly (5.16) 
\z,hc,L y ,L y = v 2L y . (5.17) 
For the width L y triangular strips, because one has to start by constructing the width- 
(L y + 1) transfer matrices, there can be more than one A even for degree d = L y . 

Concerning the A's for degree d = L y — 1, we find that one of these is the same, independent 
of L y , for the toroidal square strips we have calculated, namely 

^T,sq,L y ,L y -l,l = X . (5.18) 

This also appears for degree d = L y — 1 in the case of the cyclic square-lattice strips, and we 
conjecture that all the toroidal square-lattice strips have this eigenvalue. For the strips of 
the square lattice with L y = 3, 4, there is another common A for degree d = L y — 1, namely 

^T,sq,L y ,L y -l,2 = V ■ (5.19) 

For the toroidal honeycomb strip with L y = 4, one of the eigenvalues for degree d = L y — 1, 
say that for j — 1, is given by 

^T,hc,L y ,d=L y -l,j=l = X 2 (5.20) 

as for the corresponding cyclic strip. We conjecture that all of the toroidal strips of the 
honeycomb lattice have this eigenvalue. In terms of Potts model quantities these results are 

Xz, sg ,L y ,L y -i,i = v L y-\v + q) 

>>z,sq,L y ,L y -i,2 = v Ly {v + 1) except for L y = 2 

\ Z ,hc,L y ,L y -i,x = v 2 ^- l \v + q) 2 . (5.21) 
The expressions for the other eigenvalues are, in general, more complicated. 
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VI. SOME ILLUSTRATIVE CALCULATIONS 



A. Square-Lattice Strip, L y 



The toroidal strip of the square lattice with width L y = 2 is equivalent to the ladder graph 
with all the transverse edges doubled. The Potts model partition function Z(sq, L y x m, q, v) 
and Tutte polynomial T(sq,L y x m,x,y) were calculated for the toroidal and Klein bottle 
strips of the square lattice with width L y = 2 in Ref. 0]. We express the results here in 
terms of transfer matrices T ZjS q,2,d- For d = 0, we have 



Z,sq,2fi 



q 2 + Aqv + qv 2 + hv 2 + 2v 3 (q + 2v)(l + vf 
v 3 {2 + v) v 2 {l + v) 2 

The eigenvalues of Tz m ,i$ are the same as Xz,s2t,(5,6) given as eqs. (3.16) and (3.17) of 
For d — 1, we have 

/ 



(6.1) 



q + 3v + v 2 v(2 + v) (1+f) 2 ^ 
Tz, sq ,2,i =v v{2 + v) q + 3v + v 2 (1 + v) 2 

> v 2 (2 + v) v 2 (2 + v) v(l + vf ) 

The eigenvalues of Tz ySq ,2,i are the same as Az, S 2t,(2,3,4) given as eqs. (3.13) to (3.15) of 
The matrix TT, sq ,2,2 = v 2 has been given above. 



(6.2) 



B. Square-Lattice Strip, L y = 3 

For the L y = 3 toroidal and Klein bottle strips of the square lattice, the Potts model 
partition function Z(sq, L y xm, q, v) and Tutte polynomial T(sq, L y xm, x, y) were calculated 
in Ref. |8]. We express the results here in terms of transfer matrices T Z)Sq ,z,d- F° r d = 0, we 
have 



T: 



Z,sq,3,0 



\ 



Si VxS 2 v 1 s 2 v 1 s 2 vp 5 

V 3 S3 V 2 ViS4 V 3 ViV 2 V 3 ViV 2 V 2 v\ 

V 3 Ss V 3 ViV 2 f 2 fiS4 V 3 ViV 2 v 2 v\ 

V 3 S3 V 3 ViV 2 V 3 ViV 2 V 2 ViS4 V 2 vf 



\ V 5 V 3 V 4 V \V 2 V A ViV 2 V 4 ViV 2 v 3 v\ ) 



(6.3) 



where we use the notations v i = 1 + v , v 2 = 2 + v, v 3 = 3 + v and 



si 



q 3 + 6q 2 v + 15qv 2 



qv 



16v 3 + 3v 4 



33 



$2 = q 2 + 5qv + qv 2 + 8v 2 + 3v 3 

s 3 = q + Av + v 2 

s 4 = q + 3v + v 2 

S5 = q + 3v . 



(6.4) 



The eigenvalues of Tz tSq ,3,o consist of a linear term v 2 (v + q)(v + 1) and roots of a cubic equa- 
tion. The linear term, denoted as \z, s zt,w m H, has multiplicity two, so the corresponding 
coefficient is 26'°- ) = 2. The cubic equation is the same as eqs. (3.48) to (3.51) of [8]. For 
d — 1, we have 



T, 



Z,sq,3,l 



( 
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V 3 V iV 2 
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V 3 ViV 2 




V 3 V 1V2 





V 3 























V 2 Vi 



\ 



(6.5) 



where 



s 6 



g 2 + 5gw + 8v 2 + 



(6.6) 



The eigenvalues of Tz ySq ,3,i consist of the roots of a cubic equation, which enter with multi- 
plicity two, and the roots of a quartic equation. The cubic equation is given by eqs. (2.23) 
to (2.26) in and since it has multiplicity two, its roots have coefficients 2b^ = 2{q — 1). 
The quartic equation is the same as that given in eqs. (3.43) to (3.47) of 0. For d = 2, we 
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have 
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(6.7) 



V 

The eigenvalues of Tz, sq ^,2 consist of two linear terms and roots of three quadratic equations; 
two of the quadratic equations occur with multiplicity two. Both of the linear terms v 2 (v + q) 
and v 3 (v + 1) have the coefficient b% = (q — l)(q — 2)/2. One of the quadratic equations 
is given by eq. (3.39) of 0] with coefficient q(q — 3)/2. The other two quadratic equations 
have multiplicity two, as noted; one of them is given by eq. (2.20) of 8 ] and its roots have 



coefficient 2b, 
coefficient 2b[ 
b&. 



(2) 
2 

(2) 



{q — l)(q — 2). The other is given by eq. (2.21) of 



8| and its roots have 



q(q — 3). ThematrixTV. 



sq,3,3 



v 3 has been given above and has coefficient 



C. Square-Lattice Strip, L y = 4 



For the L y = 4 toroidal and Klein bottle strips of the square lattice, the sizes of matrices 
are 14 x 14, 35 x 35, 56 x 56, and 48 x 48 for levels < d < 3, as indicated in Table HI The 
eigenvalues of Tz, aq> A,o are roots of three quadratic equations and a sixth degree equation. One 
of the quadratic equations has multiplicity two. The eigenvalues of T Z)Sq>4:) i consist of a linear 
term v 3 (l + v)(q + v), roots of a cubic equation, a sixth-degree equation, an eighth-degree 
equation and a ninth-degree equation. The eighth-degree equation has multiplicity two. The 
eigenvalues of Tz >sq ,i } 2 consist of a linear term, roots of two cubic equations, three quartic 
equations, two sixth-degree equations and an eighth-degree equation. The linear term is 
again t> 3 (l +v)(q + v) which has multiplicity four with coefficient 2(g 2 — 3q + 1). The eighth- 
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degree equation and one of the sixth-degree equations have multiplicity two. The eigenvalues 
of T ZjS q,4:,3 consist of two linear terms t> 4 (l + v ) and v 3 (v + q), and roots of four quadratic 
equations and a quartic equation. The linear terms and one of the quadratic equations have 
multiplicity two, and the other equations have multiplicity four. The matrix Tt^q^a = v A 
has been given above with coefficient Some of these equations are too lengthy to list 
here; they are available from the authors. At the special value v — — 1, the Potts model 
partition function reduces to the chromatic polynomial. The non-zero eigenvalues are the 
same as X s tA,j for 1 < j < 33 from eqs. (2.3) to (2.21) in j^. 



D. Triangular-Lattice Strip, L y = 2 



The toroidal strip of the triangular lattice with L y = 2 is equivalent to the ladder graph 
with next-nearest-neighbor coupling and all transverse edges doubled. The corresponding 
Klein bottle strip is the same as the toroidal strip. For d = 0, we have 



T. 



q 2 + 6qv + qv 2 + 12v 2 + 8v 3 + 2v 4 (q + Av + 2v 2 )(1 + v) 



Z,tri,2,0 



v 2 (2q + I2v + I3v 2 + 6v 3 + v A ) v 2 (2 + v) 2 (l + v) 2 
The eigenvalues of Tz jiri ,2,o, denoted as Xm,{i,2)i are roots of the following equation: 



(6.S 



£ 2 - (v 6 + q 2 + 6v 5 + 20v 3 + 15w 4 + lQv 2 + qv 2 + 6vq)£ 

+v 2 (l + v) 2 {2v A + Av 3 + Av 3 q + llqv 2 + q 2 v 2 + Avq + 4q 2 v + 2q 







(6.9) 



For d = 1, we have 



/ 



Tz,tri,2,l — v 



2 \ 



q + 6v + Av 2 + v 3 q + 6v + Av 2 + v 3 (2 + v) (1 + v) 

q + 6f + Av 2 + v 3 q + Qv + Av 2 + v 3 (2 + v) (1 + v ) 2 

y v{q + I2v + 13v 2 + 6v 3 + v A ) v(q + 12v + 13v 2 + 6v 3 + v 4 ) v{2 + v) 2 {l + vf j 

(6.10) 

The eigenvalues of Tz,tri,2,i consist of a zero, Xtt2,o = 0, and roots of the following quadratic 
equation: 

£ 2 - v(v 5 + 6t» 4 + I5v 3 + 20v 2 + l&u + 2g)£ + 2v 3 (2 + v)(l + v) 2 (v 2 + vq + q) = (6.11) 
with roots A#2,(3,4)- Their coefficients are equal to = q — 1. For d — 2, we have 



rp _ 2 

1 Z,tri,2,2 — V 



1 1 
1 1 



(6.12) 
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Tz,tri,2,2 has two eigenvalues. One of them is Xtt2,5 = 2v 2 with coefficient = q(q — 3)/2; 

(2) 

while the other eigenvalue is zero with coefficient b 2 — (q — 1)(9 — 2)/2. Therefore, the Potts 
model partition function for the triangular lattice strip with L y = 2 and toroidal boundary 
condition is given by 



Z(tri,2 x L x ,tor.,q,v) = &< > X> tt2J ) L * + ^(A^-)^ + ^(^2,5) 



i=3 



(6.13) 



E. Triangular-Lattice Strip, L y = 3 

We illustrate our results for the L y = 3 toroidal strip of the triangular lattice. For d — 0, 
we have 



/ 



T; 



Z,tri,3,0 — 



V\t 2 



V\t 2 



2„,3„,2 



\ 



V t 3 V ViV 2 t4 V ViV 2 t4 V V it e V V\V 2 



V 2 V\t% V 2 ViV 2 t4 V 2 V\V 2 t4 V 2 v\v\ 



(6.14) 



where we define 



h 
U 
h 
h 



V 2 ts V 2 V\V 2 t4 V 2 V\t§ V 2 V\V 2 t4 V 2 vfv 2 

\ v A v^t-j v 3 viv 2 t 8 v 3 v\v 2 ts v 3 v\v 2 ts v 3 vfv 2 ) 



q 3 + 9q 2 v + 33qv 2 + Aqv 3 + 50v 3 + 21v 4 + 3v 5 

q 2 + 8qv + 2qv 2 + 20v 2 + Uv 3 + 3v 4 

q 2 + lOqv + 2qv 2 + 30v 2 + 22v 3 + 7v 4 + v 5 

q + Qv + Av 2 + v 3 

q + 6v + 3v 2 

2q + Uv + Uv 2 + Qv 3 + v 4 
3q + 18v + 15v 2 + 6v 3 + v 4 
q+12v + 13v 2 + 6v 3 + v 4 . 



(6.15) 



The characteristic equation of Tz,tri,z,Q yields the following cubic and quadratic equations: 

i :1 • I-mC 1 ■ hiC ■ fx; (6.16) 

where 

f 31 = -( q 3 + v 9 + 39v 2 q + 9v 8 + 36v 7 + 8Av 6 + I29v 5 + 137v 4 + 96v 3 + 12v 3 q + 2v 4 q 
+9q 2 v) 
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f32 — 



v 2 (l^ 
+2v 10 q 
+8q 3 v 5 



v)(46q 3 v 



76q 3 v 2 + 246gV + 56v 3 q 3 + 1190w 6 g + 72q 2 v 6 
v 10 + 3v u + 9v 7 q 2 + 534v 7 q + U7v 8 q + 483v 3 q 



+ q 3 v 6 



24v 9 c 



+ 9v 7 q 2 + 534v 7 q + U7v 8 q + 483v 3 q 2 + 2\3v 2 q 2 

_ 9/A. _i_ 41 A^^n _i_ _i_ 1 Q%, 4 _L 1 fil 8n 5 g _|_ 1091t> 7 



3^ "I - 35f T '■->< •<< </ II '] I hi I •/ i; "H (/ — " -i 

., , , 28gV + 1272v 4 g + 2g 4 u + 4Uv 3 q + 794v 5 + 192v 4 + 1618w 5 
+ 1243v 6 + I88v 9 + 580v 8 + 455w 4 g 2N 



-v 5 (l + vf(22q 3 v + 72gV + 105gV + 105uV + 72w 6 g + 24gV + 6g 4 + 19w 7 g 
-90u 3 g 2 + 30u 2 g 2 + 2w 4 g 4 + 12gV + 24gV + 12gV + 62gV + 60v 4 q + 19g 4 w 

4 ~ 2 ^ (6.17) 



+24w 3 g + 24w 5 + 8v 4 + 90v°q + 22v' + 30v° + 6v s + U8v*q 
with roots \ u3 j for 1 < j < 3 and 

e - v 3 (q - 2)(1 + u)£ + As - 2) 2 (1 + vf 
with roots At t3 j for j = 4, 5. For d — 1, we have 

/ 



(6.18) 
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vt u 


V 2 ViV 2 t 15 
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W Y ts 
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vt u 


V 2 *13 
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f 1)1*16 
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WiV 2 t A 


V 2 ViV 2 ti 5 


3 2 

w\v^ 


V 2 ViV 2 t 15 


v 2 t 13 




v 2 t 13 


^V^2*15 


vviv 2 U 


Wits 


WiV 2 t A 


v 2 viv 2 t 15 


3 2 

w\v 2 


VVit w 


v 2 t 12 





v 2 t l2 


f 2 f lf 2 





V 2 V\V 2 





V 2 ViV 2 








V 3 V 3 t 17 


W 3 f 3 t 17 V 3 ^17 


^1^3*18 


V 2 V 1 V 2 tg 


V 2 ViV 2 t 8 


V 2 V 1 V 2 t 8 


V 3 ViV 2 V 3 ti 8 


v 2 v\v\ 


V 3 ViV 2 V 3 t 18 





U 2 *12 


v 2 t l2 





V 2 V\V 2 








V 2 ViV 2 





V 2 ViV 2 



(6.19) 



where 



t 9 = q 2 + 8qv + qv 2 + 20v 2 + 8v 3 + v 4 

t w = 2q+ lOv + 5v 2 + v 3 

t n = q + 8v + hv 2 + v 3 

t 12 = q + 6v + 2v 2 

t 13 = 4q + qv + 24v + 20w 2 + 7v 3 + v 4 
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t u = q z + 9qv + 2qv z + 30^ + 22v 6 + 7v 4 + v 5 



t 15 = 5 + 4v + v 2 



t 16 = q + Uv + Uv 2 + 6v 3 + v 4 
t 17 = 2q + 18v + 15v 2 + 6v 3 + v 4 

t l8 = 4 + 3v + v 2 . (6.20) 



The characteristic equation of T z tri 31 yields the following quartic equation and sixth-degree 
equation: 

e 4 + hie + f^ 2 + M + fu = (6.21) 

where 

/ 41 = - v {84v 5 + 36v 6 + 9v 7 + 2v 3 q + 9v 2 q + 2q 2 + v s + UOv 3 + 130w 4 + 98v 2 + 23vq) 



/ 42 = v 3 (l + v )(Av w + 46w 9 + 2v 9 q + 2A2v 8 + 22v 8 q + 118v 7 q + 736v 7 + 383 A + 2gV 

+1388w 6 + 797v 5 q + 16gV + 1638w 5 + 1071 A + 1164v 4 + 56v 4 q 2 + 920v 3 q + 384v 3 
+109v 3 q 2 + 396v 2 q + Ulv 2 q 2 + 94q 2 v + 3q 3 v + 6q 3 ) 



/ 43 = - v Q (v + 2)(1 + v) 2 (33g 3 v + 48gV + lllgV + 42u 3 g 3 + 83Av & q + 19gV + 4v 9 q 
+3v w + 2Q2v 7 q + 4Qv 8 q + 348v 3 q 2 + 277v 2 q 2 + 3q 3 v 5 + 18q 3 v 4 + 100q 2 v + 360v 2 q 
+1836w 4 g + 1248w 3 g + 1277v 5 + 812w 4 + 12g 3 + 1585 A + 196v 3 + 602v 7 + \\27v % 
+35w 9 + 194v 8 + 268A 2 ) 

/ 44 = v 9 ( v + 2) 2 (1 + vf(Qv 6 + 18v 5 q + 10v 5 + 34v 4 q + 18v 4 q 2 + 10v 4 + 43v 3 q 2 + 4v 3 

+6v 3 q 3 + 25v 3 q + 10v 2 q + I8q 3 v 2 + 3Av 2 q 2 + I8q 3 v + lttq 2 v + 6g 3 ) (6.22) 

with roots \m,j for 6 < j < 9 and 

e + Uie + M 4 + + U£ 2 + M + fm = (6.23) 

where 

/ 61 = - v (7v 3 + 2v 4 + I2v 2 + q 2 + v 3 q + 7vq + 3v 2 q) 

/ 62 = v 2 (4t; 8 + 168w 4 + 208w 5 + l\7v & + 32v 7 + 3w 7 g + g 4 + gV + \Aq 3 v + 2w 3 g 3 + 6gV 
+73wV + 2k; 6 g + 26A 2 + 168u 4 g + 172v 3 g + hhv 3 q 2 + 6gV + 8Av 5 q) 

/ 63 = - w 6 (l + v )(2v 7 q - 8w 7 + llw 6 g + q 2 v 6 - 60v 6 + 2Av 5 q + 7q 2 v 5 - 200v 5 - 10v 4 q 
-336v 4 + 27v 4 q 2 + 2v 3 q 3 + 47w 3 g 2 - 116w 3 g - 272v 3 - 200w 2 g + 7q 3 v 2 + 32w 2 g 2 
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+I3q 3 v - 48q 2 v + q 4 - 4g 3 ) 



/ 64 = v 8 (l + w ) 2 (22g 3 w + 9gV + 38gV + 17w 3 g 3 + w 8 g 2 - 212v 6 q + 28gV + 2q 4 

+8v 7 q 2 - 6Av 7 q - 8v 8 q - 14v 3 g 2 + 104v 2 g 2 + gV + 2gV + 9gV - 4w 4 g + g 4 w 
+344w 3 g + 872v 5 + 672w 4 - 288v 5 q + 136v 7 + 492w 6 + 16v 8 + 3v 4 q 2 ) 

/ess = ^ 12 (? - 4)(u + 2)(1 + w) 3 (t; 3 + 4v 2 + Qv + q)(q + 2vf 

f m = v 14 (v + 2) 2 (1 + v) 4 (q + 2v) 4 (6.24) 

with roots X tt3 j for 10 < j < 15. Their coefficients are equal to = q — 1. For d — 2, we 
have 
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(6.25) 

where 

t 19 = q + 5v + v 2 . (6.26) 

The characteristic equation Tz,tri,z,i yields two linear equations, a quadratic equation and 
two quartic equations. The eigenvalues which are the solutions of the two linear equations, 
A«3,i6 = v 3 (v + l)(v + 2) and A tt3j i 7 = v 2 (2v + q), have the coefficient — (q — l)(g — 2)/2. 
One of the quartic equations is given by 

£ 4 - v 2 (9v 2 + 2v 3 + Uv + 2q)£ 3 

+v 4 (36v 5 + 135v 4 + 246w 3 + 192w 2 + Av 6 + 7v 3 q + 33v 2 q + 5Avq + 4g 2 )£ 2 
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-v 7 (v + 2)(1 + v)(9v 2 + 2v 3 + Uv + 2q)(q + 2v)£ + v 10 (v + 2) 2 (1 + v) 2 (g + 2v) 2 = 

(6.27) 



with roots for 18 < j < 21 that also have as their coefficient. The roots of the 

(6.28) 



following two equations have coefficient bf^ = q(q — 3)/2: 

£ 2 - v 2 (3q + 2Av + I3v 2 + 3v 3 )£ + v 5 (l + v)(6v 2 + 5vq + 6q) = 
with roots A« 3i ( 2 2,23) and 

^ + v^ 3 + v 6 (q + vq + v 2 )i 2 - v 10 q(l + v)£ + v 12 q 2 (l + vf = 
with roots for 24 < j < 27. For d = 2, we have 

( I 1 o\ 
110 



(6.29) 



Z,tri,3,3 



(6.30) 



1 1 

1 1 

1 1 
\0 1 1 / 

All the eigenvalues of T z ^h,z,z have multiplicity two. These eigenvalues are \m,2& = 2i> 3 with 
coefficient b^ = (q— l)(q 2 — 5g + 3)/3 and the roots of £ 2 — v 3 £ + v 6 = 0, denoted as A tt3i ( 2 9 i 3o), 
with coefficient 2b 2 = q{q — 2){q — 4)/3. Therefore, the Potts model partition function for 
the triangular lattice strip with L y = 3 and toroidal boundary condition is given by 

5 15 21 

Z(tri, 3 x L x , tor., q, v) = J> tt3)J ) L * + X>«3j) L " + &f E ^^^^ 



i=i i=6 

27 



i=ie 

30 



+b? E(W^ + ^ W) L * + 26® W L * . (6.31) 

j=22 j=29 

For the corresponding Klein bottle strip, the coefficients for Xm,j with 10 < j < 15, 18 < 
j < 21, 24 < j < 27 and j = 4, 5, 29, 30 become zero. The coefficient for Xtt3,n changes sign 
and the coefficient for A^ 3 2 g becomes — b^ — 1 — q. Consequently, the Potts model partition 
function for the triangular lattice strip with L y = 3 and Klein bottle boundary condition is 
given by 

3 / 9 

Z(trt,3xL x ,Kb.,q,v) = &«» ^(W^ + J> tt 3,;) L * - (Kt3,2sY 



3=1 



vj"=6 
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23 

+4 2) ((A tt 3,i6) L * - (Xm,i 7 ) Lx ) + b? ^(W 1- • ( 6 - 32 ) 

i=22 

At the special value v = —1, the Potts model partition function reduces to the chromatic 
polynomial. The non-zero eigenvalues are the same as X u j for 1 < j < 11 from eqs. (4.2) to 
(4.9) in Q 

F. Honeycomb-Lattice Strip, L y = 4 

For the L y = 4 toroidal and Klein bottle strips of the honeycomb lattice, the linear 
sizes of the matrices are 14, 35, 56, 48 for levels < d < 3, as indicated in Table The 
eigenvalues of Tz t hc,4,o are roots of two quadratic equations and an eighth-degree equation. 
One of the quadratic equations has multiplicity two. The eigenvalues of Tz,hc,±,\ are roots 
of a seventh-degree equation, an eighth-degree equation and a twelfth-degree equation. The 
eighth-degree equation has multiplicity two. The eigenvalues of Tz y hc,i,2 are roots of a quartic 
equation, a fifth-degree equation, a sixth-degree equation, two eighth-degree equations and a 
eleventh-degree equation. The sixth-degree equation and one of the eighth-degree equations 
have multiplicity two. The eigenvalues of Tzhc&s consist of a linear term v 6 (v + q) 2 , roots of 
a quadratic equation, a cubic equation and two quartic equations. The linear terms and the 
cubic equation have multiplicity two, and the other equations have multiplicity four. The 
matrix TT,he,i,& = v 8 has been given above with coefficient b^ 4 \ Some of these equations are 
too lengthy to list here. They are available from the authors. 

At the special value v — — 1, the Potts model partition function reduces to the chromatic 
polynomial. The linear sizes of the matrices reduce to 6, 18, 34, 36 for levels < d < 3, as 
indicated in Table IV II The eigenvalues of T Pjhc ^ consist of two linear terms and roots of a 
quartic equation. They are given by 

Aw4,o = (6.33) 
Am« = (?-1) 2 (?-3) 2 (6.34) 

£ 4 - (q s - 12q 7 + 66g 6 - 220g 5 + 496g 4 - 796g 3 + 922g 2 - 734g + 314)£ 3 

+ (q 12 - 20q n + 188g 10 - 1092g 9 + 4344g 8 - 12428g 7 + 26192g 6 - 41022g 5 + 47597g 4 

-40318g 3 + 24247g 2 - 9782g + 2169)£ 2 

-(4g 12 - 80g n + 736g 10 - 4124g 9 + 15705g 8 - 42922g 7 + 86535g 6 - 129964g 5 
+144575g 4 - 116374g 3 + 64525g 2 - 22280g + 3680)£ + A(q - l) 4 (q - 2) 4 = (6.35) 



42 



with roots Xhuj for 2 < j < 5. \hu,o = appears as one of the eigenvalues of Tp t h c ,4,i 
with multiplicity two. The remaining eigenvalues of T Phc ^\ are roots of a cubic equation, 
a quartic equation and a sixth-degree equation. The cubic equation has multiplicity two. 
These equations are given by 

£ 3 - (g 6 - 10g 5 + 44g 4 - 110g 3 + 167g 2 - 150g + 68)£ 2 

+ (q 8 - Uq 7 + 90g 6 - 342g 5 + 832g 4 - 1322g 3 + 1345g 2 - 812g + 230)£ 

-(q - l) 2 (2g 6 - 22g 5 + 102g 4 - 256g 3 + 373g 2 - 308g + 115) = (6.36) 

with roots \hu,j for j = 6, 7, 8. 

£ 4 - (q 6 - 8g 5 + 29g 4 - 64g 3 + 95g 2 - 86g + 37)£ 3 

+ (q - l) 2 (g 8 - 14g 7 + 86g 6 - 302g 5 + 672g 4 - lOOOg 3 + 1030g 2 - 724g + 280)£ 2 
-(? - !) 4 <y - !6g 7 + H2g 6 - 444g 5 + 1084g 4 - 1664g 3 + 1581g 2 - 886g + 253)£ 
+ ( g -3) 2 (g-l) 6 = (6.37) 

with roots Xht4,j for 9 < j < 12. 

£ 6 - (g 6 - 12g 5 + 65g 4 - 204g 3 + 407g 2 - 510g + 321)£ 5 + (g 10 - 20g 9 + 187g 8 
-1064g 7 + 4046g 6 - 10670g 5 + 19694g 4 - 25280g 3 + 22278g 2 - 13120g + 4368)£ 4 
-(g 12 - 24g n + 274g 10 - 1948g 9 + 9541g 8 - 33732g 7 + 87865g 6 - 169306g 5 + 239311g 4 
-243300g 3 + 172063g 2 - 79202g + 19177)£ 3 + (4g 12 - 92g n + 992g 10 - 6588g 9 
+29885g 8 - 97284g 7 + 232764g 6 - 412788g 5 + 540290g 4 - 512128g 3 + 336720g 2 
-139512g + 27993)£ 2 - 4(g - l) 2 (g 10 - 20g 9 + 182g 8 - 984g 7 + 3483g 6 - 8402g 5 
+ 13956g 4 - 15756g 3 + 11602g 2 - 5078g + 1032)£ + 16(g - 2) 2 (g - l) 6 = (6.38) 

with roots Xht4,j for 13 < j < 18. Xht4,o = is again one of the eigenvalues of Tp^ c ,4,2- The 
remaining eigenvalues of T Pjhc ^ 2 are roots of a quadratic equation, two cubic equations, two 
fifth-degree equations and a seventh-degree equation. They are given by 

A W 4,(i9,20) = \ [q 2 - 4g + 7 ± (g - 3)(g 2 - 2g + 5) 1 / 2 ] (6.39) 

£ 5 - (2g 4 - 14g 3 + 41g 2 - 60g + 41)£ 4 + (g 8 - 14g 7 + 91g 6 - 358g 5 + 932g 4 - 1636g 3 
+ 1875g 2 - 1270g + 396)£ 3 - (g 10 - 18g 9 + 148g 8 - 726g 7 + 2339g 6 - 5148g 5 + 7824g 4 
-8136g 3 + 5615g 2 - 2372g + 477)£ 2 + (g - l) 2 (g 8 - 14g 7 + 85g 6 - 292g 5 + 626g 4 
-874g 3 + 800g 2 - 448g + 125)£ - 4(g - l) 4 = (6.40) 

with roots X hU j for 21 < j < 25. 

£ 5 - (2g 4 - 15g 3 + 51g 2 - 90g + 73)£ 4 + (g 8 - 15g 7 + 104g 6 - 429g 5 + 1152g 4 - 2064g 3 
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+2415g 2 - 1694g + 570)£ 3 - (g 10 - 19g 9 + 167g 8 - 881g 7 + 3070g 6 - 7372g 5 + 12396g 4 
-14520g 3 + 11459g 2 - 5550g + 1265)£ 2 + (g - l) 2 (g 8 - 15g 7 + 104g 6 - 428g 5 + 1144g 4 
-2039g 3 + 2392g 2 - 1706g + 575)£ - 2g(g - l) 6 = (6.41) 

with roots Xht4,j for 26 < j < 30. 

£ 3 - (q 4 - 7q 3 + 21g 2 - 32g + 24)£ 2 + (g 6 - llg 5 + 53g 4 - 139g 3 + 206g 2 - 164g + 58)£ 
-(? - l ?{(f ~ 7q 3 + 18g 2 - 24g + 15) = (6.42) 

with roots Xhuj for j = 31, 32, 33. 

£ 3 - (g 4 - 8g 3 + 27g 2 - 40g + 24)£ 2 + (g - l) 2 (g 4 - 10g 3 + 44g 2 - 90g + 72)£ 
-(g-3) 2 (g-l) 4 = (6.43) 

with roots \hu,j for j = 34, 35, 36. 

£ 7 - (3g 4 - 22g 3 + 75g 2 - 140g + 126)£ 6 + (3g 8 - 44g 7 + 304g 6 - 1282g 5 + 3601g 4 
-6890g 3 + 8830g 2 - 7038g + 2781)£ 5 - (g 12 - 22g n + 230g 10 - 1504g 9 + 6839g 8 
-22756g 7 + 56785g 6 - 107160g 5 + 152314g 4 - 160280g 3 + 120302g 2 - 59250g + 15045)£ 4 
+ (g 14 - 26g 13 + 318g 12 - 2414g n + 12685g 10 - 48806g 9 + 141948g 8 - 317702g 7 
+551665g 6 - 743122g 5 + 769969g 4 - 601782g 3 + 341290g 2 - 128886g + 25246)£ 3 
-(g 14 - 24g 13 + 280g 12 - 2088g n + 11028g 10 - 43286g 9 + 129310g 8 - 297078g 7 
+525432g 6 - 710236g 5 + 721717g 4 - 535598g 3 + 276095g 2 - 89634g + 14145)£ 2 
+4(g - l) 2 (g 9 - 6g 8 - 14g 7 + 249g 6 - 1080g 5 + 2522g 4 - 3564g 3 + 3069g 2 - 1499g 
+326)£-4(g-2) 2 (g- l) 6 = (6.44) 

with roots \hu,j for 37 < j < 43. The second fifth- degree equation and the first cubic 
equation given above have multiplicity two. A^o = is again one of the eigenvalues of 
Tp,hc,4,3 with multiplicity two. The remaining eigenvalues of Tp^c,^ are two linear terms, 
roots of a quadratic equation and two cubic equations. They are given by 

A hi4i 44 = g 2 - 4g + 5 (6.45) 

A ht4) 45 = (q - I) 2 (6-46) 

A W4 ,(46,47) = \ [q 2 ~ 6g + 13 ± (g - 3)(g 2 - 6g + 17) 1 / 2 ] (6.47) 

£ 3 - (2g 2 - 8g + 11)£ 2 + (g 4 - 8g 3 + 26g 2 - 38g + 23)£ - (g - l) 2 = (6.48) 
with roots \hu,j for j = 48, 49, 50. 

£ 3 - (2g 2 - 8g + 13)£ 2 + (g 4 - 8g 3 + 26g 2 - 34g + 19)£ - 3(g - l) 2 = (6.49) 
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with roots \ht4,j for j = 51, 52, 53. Xht4,j for j = 45, 46, 47 have multiplicity two, and the other 
eigenvalues of T Phc ^ 3 have multiplicity four. Finally, we have Xht4,54 = 1 as the eigenvalue 
of Tp t hc,A,4- The corresponding coefficients are 

ht4,j = 1 for 1 < j < 5 

b ht 4j = 2{q - 1) for 6 < j < 8 

6ftt4j =9-1 for 9 < j < 18 

&M4j = ^(9 - - 2) for 19 < j < 25 

6fct4j = (? - - 2) for 26 < j < 30 

^*4,i = <?(<? - 3) for 31 < j < 33 

&«4,i = \q{q - 3) for 34 < j < 43 
2 

&M4.44 = 1)(? 2 - 5g + 3) 

6w4j = - 1)(? 2 - 5g + 3) for 45 < j < 47 
&ftt4j = \q{q ~ 2)(q - 4) for 48 < j < 53 

b h u M = g 4 - 8g 3 + 20g 2 - 15q + 1 . (6.50) 

Although XhtAfl = does not contribute the the chromatic polynomial, we can still calculate 
its coefficient for the toroidal strip, and we get (q — l)(2g 2 — 7q + 12)/6. The sum of all the 
coefficients is equal to q 2 (q — l) 2 , as dictated by the L y = 4 special case of the general result 
(|3.14|) . Hence, the chromatic polynomial for the honeycomb lattice strip with L y = 4 and 
toroidal boundary condition is given by 

5 / 8 18 

P(hc,4xL x ,tor.,q) = b® ^(W™ + 2 ^(^W m + ^(A W4j ) 

3=1 V J=6 i=9 

( 25 30 
^(A, t4j r + 2^(A, t4 , J ) 
j=19 i=26 



45 



33 43 

m 



-bf^ J 2 ^^(Aht4j) m + ^(\(4j, 

V i=31 i=34 
/ 47 

-fe^ J 2(A/ lt 4 > 4 4 ) m + y^(Afef4j) m 

j=45 

53 



+44 3) ^(A, t4j ) m + 6 (4) (AM4,5 4 ) m (6.51) 

i=48 

where m = L x /2 as given in eq. (jl.4j) . For the corresponding Klein bottle strip, Xhu,j for 
6 < j < 8, 26 < j < 33, 48 < j < 53 and j = 44 do not contribute. The coefficient for 
Am4j with 9 < j < 12, 21 < j < 25, 34 < j < 36 and j = 1 changes sign. The coefficient 
for Am4,45 becomes = q — 1, the coefficients for \hu,[m,A7) become — feW — 1 — q and the 
coefficients for A^4,54 become one. Therefore, the chromatic polynomial for the honeycomb 
lattice strip with L y = 4 and Klein bottle boundary condition is given by 

P(hc,AxL x ,Kb.,q) = &<°> (-(A M4 ,i) m + ^(A^) m + (Am4, 54 ) 

V 3=2 

12 18 47 

- ^2(^hu,j) m + y^(AM4j) m + (A/i t 4,45) m - y^(A/ t M,j) m 

i=9 j=13 i=46 

20 25 



,m 
,3! 



36 43 



\ j=34 i=37 / 

The singular locus B for the L x — > oo limit of the strip of the honeycomb lattice with 
L y = 4 and toroidal boundary conditions is shown in Fig. For comparison, chromatic 
zeros are calculated and shown for length L x = 13 (i.e., n = 104 vertices). The locus B 
crosses the real axis at the points q = 0, q = 2, and at the maximal point q = q c , where 

q c ~ 2.480749 for {G} = (he, 4 x oo, PBC y , PBC X ) . (6.53) 

We note that this is rather close to the value for the infinite 2D honeycomb lattice, namely 
q c (hc) = (3 + \/5 )/2 ~ 2.61803. As is evident from Fig. [TUl the locus B separates the q 
plane into different regions including the following: (i) region Ri, containing the semi-infinite 
intervals q > q c and q < on the real axis and extending outward to infinite \q\; (ii) R2 
containing the interval 2 < q < q c ; (iii) i? 3 containing the real interval < q < 2; and (iv) 
the complex-conjugate pair i? 4 , R* A centered approximately at q = 2.55 ± 0.6i. The (nonzero) 
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FIG. 10: Singular locus B for the L x — > oo limit of the strip of the honeycomb lattice with L y = 4 and 
toroidal boundary conditions. For comparison, chromatic zeros are shown for L x = 13 (i.e., n = 104). 

densities of chromatic zeros have the smallest values on the curve separating regions R± and 
i?3 in the vicinity of the point q = and on the curve separating regions R2 and R3 in the 
vicinity of the point q = 2. 

It is instructive to compare this locus B in Fig. 1101 with the corresponding loci that we 
obtained previously for strips of the honeycomb lattice with periodic longitudinal boundary 
conditions and free or periodic transverse boundary conditions. We recall that for a given 
infinite-length lattice strip with specified transverse boundary conditions, the locus B is 
the same for periodic longitudinal and twisted periodic longitudinal boundary conditions 
(i.e., cyclic versus Mobius if the transverse B.C. are free, and toroidal versus Klein bottle 
if the transverse B.C. are periodic). For the cyclic strip of the honeycomb lattice, B was 



given as (i) Fig. 1 of Ref. |25j for width L y = 2± and (ii) Figs. 7, 8, and 9 of |22j for 
L y = 3,4,5 (see also the zeros in Fig. 17 of Ref. |33J] for L y — 3). We recall that L y must 
be even for strips of the honeycomb lattice with toroidal boundary conditions. Comparing 
our results for B for the width L y = 4 toroidal strip of the honeycomb lattice in Fig. El 



221, we see 



with our earlier results for the L y = 4 cyclic strip of this lattice in Fig. 8 of 
that the locus B is significantly simpler for toroidal, versus cyclic, boundary conditions, with 
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fewer regions and, in particular, no evidence of the tiny sliver phases that are present in the 
latter case. Furthermore, the value of q c ~ 2.481 that we obtain (eq. (|6.53Jl ) for the width 
L y = 4 infinite-length strip of the honeycomb lattice with toroidal boundary conditions is 
closer to the value of q c {hc) ~ 2.618 for the infinite 2D than the value q c ~ 2.155 for the 
corresponding cyclic L y = 4 strip. This is in agreement with one's general expectations, 
since toroidal B.C. minimize finite-size effects to a greater extent than cyclic B.C. and hence 
expedite the approach to the L y — ► oo limit. 

In region R\, the dominant eigenvalue is the root of eq. (|6.35j) with largest magnitude. 
Denoting it as Xhu,2, we have 

W = (X ht4 , 2 ) 1/8 , qeRt. (6.54) 

This is the same as W for the corresponding L x — > oo limit of the strip of the honeycomb 
lattice with the same width L y = 4 and cylindrical (PBC y , FBC X ) boundary conditions, 
calculated in [3^. This equality of the W functions for the L x — > 00 limit of two strips 
of a given lattice with the same transverse boundary conditions and different longitudinal 
boundary conditions in the more restrictive region Ri defined by the two boundary conditions 
is a general result 0, |45|. In we have listed values of W for a range of values of q for the 
L x — > 00 limit of various strips of the honeycomb lattice, including (he, 4x 00, PBC y , FBC X ). 
Since W is independent of BC X for q in the more restrictive region Ri defined by FBC X and 
(T)PBC X (which is the R\ defined by PBC X here), it follows, in particular, that 

W(A x 00, PBC y , (T)PBC X , q) = W(A x 00, PBC y , FBC X , q) for q > 2.48... (6.55) 

In region R 2 , the largest root of the seventh degree equation (|6.44j) is dominant; we label 
this as \ha,37, so that 

\W\ = \X h u,3 7 \ 1/8 , q e R 2 (6.56) 
(in regions other than R x , only \W\ can be determined unambiguously [7]). Thus, q c is the 
relevant solution of the equation of degeneracy in magnitude (A^^l = ^44,37 1- In region 
R3, 

\W\ = |A ht4il3 r /8 , q e R 3 (6.57) 
where Xhu,i3 is the largest root of the sixth degree equation (|6.38|) . In regions _R 4 , RI, 

\W\ = \X ht4 , 9 \ 1/8 , qeRt,Rt (6.58) 

where g is the largest root of the quartic equation (|6.37|) . We calculate the following 
specific values: \W{q = 1)| = 1.67290.., \W(q = 2)| = 1.35817.., and \W(q = q c )\ = 1.26949.. 
These values may be compared with those listed in Ref . [3^ . As is evident from Fig. the 
curve B has support for Re(q) < 0. For strips with finite lengths m > 3, certain chromatic 
zeros have support for Re(q) < 0. 
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VII. SUMMARY 



In summary, in this paper we have presented a method for calculating transfer matrices 
for the g-state Potts model partition functions Z(G,q,v), for arbitrary q and temperature 
variable v, on strip graphs G of the square, triangular, and honeycomb lattices with width 
L y and arbitrarily great length L x , having toroidal and Klein bottle boundary conditions. 
Using this method, we have derived a number of general properties of these transfer matrices. 
In particular, we have found some very simple formulas for the determinant det{Tz,K,L y ,d), 
and trace Tr(T Z A Lj/ ). A number of explicit exact calculations were given. 
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